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DEFINITION OF SYMBOLS 


Symbol 

A 0 

B N 

B P 

b 0 

b NP 

C N-1,N 

F 

k n 

k 

N 

n 

P 

R 

S NN 

S NY 

S YY 

S(Reg) 

S(Res) 

s 

Y 


Definition 

model equation intercept coefficient 

matrix quantity 

model equation coefficients 

defined quantity 

regression model coefficients 

regression model cross product coefficients 

F statistic 

cross product model equation terms 

cross product technique parameter pertaining to degrees of freedom 

number of independent variables 

number of data points 

degree of equation 

multiple correlation coefficient 

matrix quantity 

matrix quantity 

total sum of squares 

regression sum of squares 

residual sum of squares 

variance 

value of dependent variable 


v 



DEFINITION OF SYMBOLS (Concluded) 


Symbol 

y 

Z N 

Z NP 


AVG 



Definition 

difference in observed and mean value of dependent variable 
value of Nth independent variable 

difference in observed and mean values to the Pth power for the Nth 
independent variable 

average absolute percent error 

n 



, summation from the first to the nth term 


Superscript 

c calculated value 

— mean value 


Subscript 

i i = 1 , 2, ... n 

j j = 1 , 2, ... N 


vi 



MULTIPLE REGRESSION TECHNIQUE FOR Pth DEGREE POLYNOMIALS 
WITH AND WITHOUT LINEAR CROSS PRODUCTS 

l. INTRODUCTION 


The investigation of physical processes and requirements for data analysis methods 
frequently requires the use of mathematical models which describe the processes. The 
model can be formulated such that certain variables interact according to physical 
theories associated with the particular process, or it may contain identified independent 
variables and unknown parameters. The relationship of the parameters identified in the 
model can be evaluated using the statistical tool commonly referred to as regression 
analysis. In principle, it should be possible to establish complex curves or surfaces for 
higher order multiple variable functions by regression techniques to summarize trends in 
data and to provide a means of predicting similar phenomena. Furthermore, such a 
technique might be used to establish unknown laws or relationships. 

Most statistical textbooks treat the problem of linear multiple variable regression 
and of nonlinear regression of one independent variable. However, the study of many 
complex physical problems requires a method capable of determining nonlinear regression 
of multiple independent variables. By this means, an analytical representation of the 
observed data is provided which can be used in subsequent analysis. 

The general procedure in regression analysis is to take partial derivatives of a 
specific model-dependent minimizing function. The set of equations obtained by setting 
these partial derivatives equal to zero is frequently referred to as the normal equations. 
If the normal equations are not transcendental in any of the unknown parameters, they 
can be solved by the usual algebraic methods. It is this situation which is of concern 
here. 


Within this report a method is presented for evaluating multiple variable regression 
for Pth degree polynomials with and without linear cross products. Implementation of 
the technique, including error diagnostics has been accomplished on the UNIVAC 1108 
computer. Program listings and illustrative example problems are given in the appendices. 

The method developed is used in problems requiring the determination of 
mathematical relationslups describing complex curves or surfaces from known physical 
data. Initial applications of this technique have concerned the external aerodynamics of 
aircraft and space vehicles and the internal aerodynamics of transonic wind tunnels and 
have produced excellent results. Applications in almost any other field of study should be 
expected to be equally as useful. 



II. MULTIPLE REGRESSION TECHNIQUE WITHOUT CROSS PRODUCTS 


A. Introduction 

The derivation of a matrix solution for evaluating multiple variable regression for 
Pth degree polynomials without cross products is presented in Reference 1. In this 
section this solution is discussed for the general case where the polynomial expression can 
be of Pth degree with N independent variables. The mathematical procedures represent a 
rigorous least squares evolution of input data centered about the mean. 


8. Regression Analysis Development 

Assuming the observed dependent variable is to be estimated by the following 

model: 


Y i A 0 + ^ll Z li + b 12 Z li + ‘ ' + b,D Z,: ^ 


IP li 


+ (b 2l Z 2i + b 22 Z 2i + -‘ +b - Z m ) 


2P 2i 


+ . . . + (b N1 Z Nj + b N2 ZJjj + . . . + b Np Z N . ) 


( 1 ) 


where i = 1 , 2, nth set of data, and the intercept = b 1 A + b ™ + ... + b XT ~ . 

r 0 10 20 NO 

For any given independent variable , the mean value is given by 


z 

N n 

It the input data are centered about the mean the model equation (1) becomes 
Y i = b 0 + (b l 1 Z lli + b 12 z 12i + + b lP z lPi ) 

+ ( ' b 21 z 21i + b 22 Z 22i + ’ ' + b 2P z 2Pi * 


( 2 ) 


(3) 


2 




Where, for convenience, the following quantities, which are indicated in equation (3), are 
defined by 


b 0 (b 10 + b l 1 Z 1 + b 12 V + ' ' ' + b lP Z ! } 


+ (b 2Q + b. 21 Z 2 + b 22 Z 2 2 + . . . + b 2p Z 2 ) 



(4) 
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The classical form of the least squares minimizing function is 

M - S CYj - Yj 0 ) 3 (6) 


This result follows from the least squares principle that the best representation of the 
input data is that which makes the sum of the squares of the residuals a minimum. The 
condition which fulfills this requirement is that the partial derivatives of this function 
with respect to each of the unknowns be zero. Hence, the following normal equations 
for i = 1 , 2, ... n are written: 


3M 

3b„ 


= 0 


3M 


8b 


- 0 


1 1 


8M 
3 1 2 


- 0 


3M 


3b 


= 0 


IP 


9M 
3 b 


21 


3M 


3 b 


= 0 


- 0 


22 


3M 


3b 


= 0 


2P 


3M 


3 b 


N1 


3M 


3b 


N2 


3M 


3b 


= 0 


= 0 


= 0 


NP 


(7) 


Through the use of the normal equation ~j~ = 0 it can be shown that 

3b 0 


b 0 ~ Y = 




m 


( 8 ) 


The following quantity is now defined for the difference in the observed value 
of Yj and the mean value of the observed values: 


Y. - Y 


(9) 
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Then, through the use of the remaining normal equations, the following matrix 
solution for the unknown constant of the independent variables is obtained. 



where the following matrix quantities are given by: 



5 



’lj 


2j 


Z z iii z jii 

Z Z Hi Z j2i - 

Z z l li Z jPi 

Z Z 12i Z jli 

Z Z 12i Z j2i • • 

Z Z 12i z jPi 

Z Z 1 Pi Z jli 

Z Z lPi Z j2i 

Z Z lPi Z jPi 

Z Z 2 1 i Z j 1 i 

Z Z 21i Z j2i 

Z Z 21i Z jPi 

Z z 22i Z jli 

Z Z 22i Z j2i 

Z z 22i Z jPi 

Z Z 2Pi z jli 

Z Z 2Pi Z j2i 

Z Z 2Pi z jPi 



for j = 1, 2, N 


for j = 1, 2, N 





Z Z Nli Z jli 
S Z N2i z jli 

Z Z NPi z jli 


E Z Nli z j2i 

Z Z N2i Z j2i 

Z z NPi z j2i 


E Z Nli Z jPi 
E Z N2i Z jPi 

Z Z NPi Z jPi 


for j = 1, 2, N 

J 


( 12 ) 
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Having determined the matrix solution indicated by equation (10), the intercept 

A of the fitted model expression may be determined as: 
o 

A 0 = Y_ (b ll Z 1 +b 12 Z \ + - • • + b ip^i )_(b 21 Z 2 + h 22 z 2 + ‘ " 

+ b 2 Z 2 P)- '- “ (b Nl Z N + b N2 Z N + -- < + b NP Z N ) (14) 
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C. Discussion 


As noted by Graybill [2], there are an infinite number of solutions to the model 
relationship. However, only one solution must be found to have a useful result. Also note 
that the solution presented to the higher order multiple regression problem is a function 
of parameters such as the sum of the squares, cross products, and linear values similar to 
previously developed solutions for linear regression. In this case, one finds a more 
complex result with nested matrices, which becomes practical only in combination with 
digital computer techniques. 

As discussed in Appendix A, a computer program has been developed to evaluate 
the unknown constants in the model equation with the solution specified by equation 
(10), as well as certain diagnostics reflecting the significance of the estimated regression 
equation which will be considered in Section IV. 

It is required in the application of this technique that the S matrix be nonsingular 
and that the input data are reasonably well behaved. It is also required that n = NP . 
That is, the number of data points n must be equal to or greater than the number of 
unknown parameters in the model equation. 

Many sets of physical data have been studied using the previously mentioned 
computer program. Results have generally been excellent and it is felt that the technique 
can be highly useful in many fields of endeavor. 


III. MULTIPLE REGRESSION TECHNIQUE WITH CROSS PRODUCTS 


A. Introduction 

The regression technique developed in Section II for Pth degree polynomials with 
N independent variables has proven to be extremely useful in correlating experimental 
wind tunnel data and aiding subsequent analysis. However, in certain applications the 
inclusion of cross products for the various independent variables can enhance the value of 
the technique. Such a case occurs when it is desired to optimize the dependent variable 
of the fitted mathematical expression with respect to the various independent variables. 
In such a case the determination of maximum and minimum values by solving for the 
critical points through the use of the necessary and sufficient conditions (that the first 
and second partial derivation of the response variable with respect to each of the 
independent variables be identically equal to zero) can only be determined when cross 
products are included in the model equation; otherwise the effect of all other 
independent variables is lost when the partial derivative with respect to a given 
independent variable is evaluated. 
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To deal with this and other similar problems, a method is developed in this 
section which treats all combinations of linear cross products for each of the N 
independent variables as well as the polynomial terms previously discussed in Section II. 
Higher order cross product terms could, in principle, be treated, although such a 
development is beyond the scope of the present work. The method for fitting 
polynomials with linear cross products is presented without formal derivation. However, 
the development is similar to that shown in Reference 1. 


B. Regression Analysis Development with Linear Cross Products 

Consider the following model: 


n 

> 

o 

+ 

(bll 

z ii 

+ b 12 Z li + ' 

. , + bjp 

Z ^) 
4 li ] 

+ 

(b 21 

Z 2i 

+ b 22 Z 2i + ' 

. . + b 2p 

Z ^) 
^2i ] 

+ . . . + 

< b Nl 

Z Ni 

+ b N2 Z Ni + 

••• + b NP Z Ni P) 

+ 

(C 12 

Z li 

Z 2i + C 13 Z 1 

. Z-. . 

i 3i 

. + c lN Z i 

+ 

(C 23 

Z 2, 

Z 3i + C 24 Z 

Z A . + . 
2i 4i 

. . + c 2N Z, 


+ -' +(C N-l,N Z (N-l)i Z Ni } 


(15) 


where i = 1 , 2, ... , nth set of data, and the intercept Aq = b j q + b 2 Q + ... + b n g . This 

model represents a polynomial expression of Pth degree with N independent variables and 
contains all combinations of linear cross products for each of the N independent 
variables. 

If equal quantities are added and subtracted in equation (15), namely, those of 
the form bj^jpZ^ and CjjZjZj , (where j = 1, 2, N-l) it can be shown that the model 
equation is equivalent to: 
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(16) 


Y; = K, • + B , . + K-. + B«. + . . . + K XT . + B x 


L . lVi . < LJ i ■ 1 lV,v 1 LI . 

i h li 2i 2i 


Ni Ni 


where 


K li b 10 + b ll Z 1 +b 12 Z l + • +b lP Z l +C 12 Z I Z 2 +C 13 Z l Z 3 + " +C lN Z l Z N 


— 2 


~ P 


K 2i b 20 + b 21 Z 2 + b 22 Z 2 + ,, ‘ + b 2P Z 2 +C 23 Z 2 Z 3 +C 24 Z 2 Z 4 + ' ' ’ +C 2N Z 2 Z N 1 


— 2 


7 P 


K xt , - b XTn + b XT1 Z XT + b XT ^Z VT + . .. + b Np Z N +C n_i 


k Ni NO NI N N2 N 


(17) 


B li b ll Z lli + b 12 Z 12i + - +b lP Z lPi +C 12 Z 12i +C l3 Z 13i + - ’ +C lN Z 1 Ni 


B 2i b 21 z 21i + b 22 z 22i + ' ' ' + b 2P Z 2Pi +C 23 Z 23i +C 24 z 24i + ’" +C 2N Z 2Ni 


B Ni b N 1 z Nli + b N2 Z N2i + • + b NP Z NP +C N-1 ,N Z N-l,Nj 


(18) 


The jth mean value terms in equation (18) are given by 


z. = T JL 

L ) L* n 


(19) 




( 20 ) 


and for convenience z and z terms in equation (18) have been defined as: 
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z lli ' z li‘ z l 
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Z 23i Z 2i Z 3i Z 2 Z 3 


z t A . = Z~. Z.. - ZJL . 
24i 2i 4i 2 4 


z+ = Z 0 . Z M . - Z_Z. T 
2Ni 2i Ni 2 N 


> 


( 22 ) 

(Concluded) 



7 7-7 7 

N-l,i Ni N-1N 



Regrouping equation (2) the following is obtained 



b 0 + B ii + B 2i + 


B 


Ni 


(23) 


where 


b 0- K U + K 2i + - + K Ni 


(24) 


Through the use of the normal equations resulting from the least squares 
minimizing function it is found that 



(25) 


The difference in the observed value and the mean value of the observed values of 
the dependent variable is again denoted by: 


y i = 


Y. - Y 
x 


(26) 
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Then the following matrix .solution for the unknown constants of the model 
equation can be obtained through the use of the normal equations: 


r h 


— 



-1 

p — 

B i 


S 11 

S 12 • 

S 1N 


S 1Y 

B 2 

= 

S 21 

S 22 ' ' 

S 2N 


S 2Y 

1 

55 

z • • 

1 


S N1 

S N2 • 1 

• S NN 


S NY 


[Mxl] [MxM] [Mxl] 

where 

M = NP + N!/2(N— 2)! 


(27) 


(28) 


and the various matrix quantities are given by: 


B, = 


B„ = 


" b iT 


b ip 


C 12 


C 1N 


" b 2 r 


b 2P 


c 23 


_ C 2N_ 





(29) 
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14 



S 12 " 


S 22 " 


S N2 “ 


E z lli z 21i ••• 

E Z 1 li z 2Pi 

Z Z 1 Pi z 21i ••• 

E z lPi z 2Pi 

Z z 12i z 21i ••• 

E z 12i z 1 Pi 

Z z *Ni z 21i 

E z l*Ni z 2Pi 

Z z 2Ii z 21i ••• 

E z 21i z 2Pi 

Z z 2Pi z 21i ••• 

E z 2Pi z 2Pi 

Z z 23i z 21i ••• 

E z 23i z 2Pi 

Z z 2Ni z 21i ••• 

« 

• 

E z 2Ni z 2Pi 

♦ 

Z z Nli z 21i 

E z Nli z 2Pi 

Z z NPi z 21i ••• 

E z NPi z 2Pi 

Z z N-l,Ni z 21i • 

• • E Z N-1 ,Ni z 2Pi 


E z lH z 23i 

E Z 1 li z 2Ni 

Z z lPi z 23i 

E z 1 Pi z 2Ni 

Z z 12i z 23i ••• 

E z 12i z 2Ni 

Z z lNi z 23i ••• 

z lNi z 2Ni 

E z 21i z 23i 

E z 21i z 2Ni 

E z 2Pi z 23i • • ■ 

E z 2Pi z 2Ni 

E z 23i z 23i • • • 

E z 23i z 2Ni 

E z 2Ni z 23i • • • 

ri + + 

z 2Ni z 2Ni 

E z Nli z 23i ■ ■ • 

E Z N 1 i z 2Ni 

E z NPi z 23i ■ • • 

E z NPi z 2Ni 

E Z N-1 ,Ni z 23i • • 

v + _+ 

■ • h z N-l,Ni z 2Ni 







(30) 

(Continued) 
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*\ 


Z z lli z Nli ••• 

Z Z 1 li z NPi 

Z Z 1 li z N-l,Ni 

' Z Z 1 Pi z Nli • • • 

Z z 1 Pi z NPi 

Z z lPi z N-l,Ni 

Z z i2i z Nii ••• 

Z z 12i z NPi 

Z z ^2i z N-l,Ni 

Z z lNi z Nli • • • 

Z z lNi z NPi 

Z z INi z N-l,Ni 

Z z 21 z Nli 

Z z 21i z NPi 

Z z 21i z N-l,Ni 

Z z 2Pi z Nli • • • 

Z z 2Pi z NPi 

Z z 2Pi Z N-1 ,Ni 

Z z 23i z Nli • * • 

Z z 23i z NPi 

Z z 23i z N-l,Ni 

Z z 2Ni z Nli • ■ • 

• 

« 

Z z 2Ni z NPi 

h z 2Ni z N-l,Ni 

• 

Z Z N 1 i Z N 1 i 

Z z Nli z NPi 

Z Z N 1 i z N-l,Ni 

Z z NPi Z N 1 i • • • 

Z z NPi z NPi 

Z z NPi z N-l,Ni 

Z z N-l,Ni z Nli • 

• 'Z z N-l,Ni z NPi 

Z Z N-1 ,Ni Z N-1 ,Ni 




(30) 

(Concluded) 


16 




Z z 21i y i 

Z z 2Pi y i 

Z z 23i y i 
Z z 2Ni y i 


Ny 


Z z Nli y i 
Z Z NPi y i 









(31) 


Once the matrix solution given by equation (27) has been determined the 
intercept Aq of the fitted model expression may be calculated as 
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A 0 = Y-(b u Zj +b 12 Zj 2 + ... + b Jp Z 1 +C 1 2 Z 1 Z 2 + C 13 Z 1 Z 3 + ••• + c in Z 1 Z n) 


- (b 21 Z, + b 22 Z 2 + . . . + b 2p Z 2 + c 23 Z 2 Z 3 + C 24 Z 2 Z 4 + • • • + c 2N Z 2 Z n * 


~ P 


(b Nl +b N2 Z N + • * ' + b NP Z N +C N-1,N Z N-1 Z N } 


( 32 ) 


C. Comments 

The solution developed for evaluating multiple variable regression for Pth degree 
polynomials with linear cross products is similar to that developed in Section II for the 
case without cross products. However, in this case, the nested matrices contain additional 

> N! 

terms to account for the cross products and hence it is required that n = NP + ^ , 

which again indicates that the number of data points n must be equal to or greater than 
the number of unknown parameters in the model equation. 

As discussed hi Appendix B, a computer program has been developed to treat this 
case. It should be noted that while contrived problems have yielded excellent results, 
experience with several sets of physical data have generally yielded poorer results. This is 
evidently due to the larger matrices being manipulated and due to unknowns as to 
whether true linear cross coupling exists in the given physical process. 


IV. SIGNIFICANCE OF THE ESTIMATED REGRESSION EQUATIONS 


A, Introduction 

To determine if the fitted regression equation obtained from the solution matrix 
is indeed a useful representation of the input data, it is desirable that certain diagnostics 
be evaluated. To this end the computer programs discussed in Appendix A and Appendix 
B determine the standard deviation of the observed data with respect to the fitted 
equation, the multiple correlation coefficient, and the F statistic as well as the average error 
and the maximum error of the observed data with respect to the fitted result. 

Using these parameters, it is possible to assess the usefulness of the fitted expression 
for each given application. 
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B. Analysis of Variance 


The significance of the estimated regression equation can be considered from the 
viewpoint of an analysis of variance as summarized in Table 1, where the total sum of 
squares is resolved into a component measuring the residual fitting error, and a 
component which measures the regression variation being tested. 


TABLE 1. ANALYSIS OF VARIANCE (WITH AND WITHOUT CROSS PRODUCTS) 
A. Analysis of Variance (without Cross Products) 


Degrees of 
Freedom 

Type 

Variation 

Sum of Squares 
(SS) 

Mean Square 
(MS) 

F 

Value 

n-l 

NP-1 

n-NP 

Total 

Residual 

Regression 

S = Z(Y.-Y) 2 

yy 1 

S(RES) = £(Y.-Y. C ) 2 

S(REG) = ^(Y. C - Y) 2 

S(RES) 

M(RES) “ NP-l 

S(REG) 

M(REG) 

M(KhCj) 

n-NP 

M(RES) 


B. Analysis of Variance (with Cross Products) 


Degrees of 
Freedom 

Type 

Variation 

Sum of Squares 
(SS) 

Mean Square 
(MS) 

F 

Value 


Total 

s yy = Z<Vj-Y) ! 



K-l 

Residual 

S(RES) = E (Y i" Y i C)2 

S(REG) 

M(RES) = 

k-l 


n-k 

Regression 

S(REG) = ^(Y^-Y) 2 

S(RES) 
M(REG) = — ■ — 
n-k 

M(REG) 

M(RES) 


where MS = SS/degrees of freedom 


S = S(RES) + S(REG) 

yy 

Y — average of observed values 
k = NP + (N-l ) + (N-2) + . . . + 1 
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C. Interpretation Diagnostics 


As noted by Smille [3], the results of such an analysis of variance can be used to 
test the combined effect of all of the independent variables on the dependent variable. 
That is, the hypothesis that all of the population regression coefficients in the model 
regression equation are zero can be tested since the ratios ot the regression mean square 
to the residual mean square are distributed in an F distribution as shown below: 


S(REG)/Regression degrees of freedom 
S(RES)/Residual degrees of freedom 


(33) 


where it is assumed that the observations are selected at random- from a normally 

distributed population with zero mean and constant variance, that Zj , , ... are 

independent variables following x 2 distributions, and that only random errors are 
associated with the observations. 

The F ratio calculated from equation (33) can be used to test the statistical 

significance of the regression equation under consideration by comparing it with the 

appropriate FTABLE va ^ ue desired probability level with the specified numerator 

and denominator degrees of freedom; that is, the following test of the null hypothesis 
may be performed: 


1) 

o 

xT 

o 

X 

b 12 = ••• b NP = 0 


accept when 

F c > F TABLE 

(34) 

reject when 

F c < F TABLE 



D. Output Data 

The desired coefficients of the model equation are outputted in E notation 
where Ao is the computed intercept of the fitted polynomial expression and the B 
coefficients are printed out in ascending order of degree with the first P coefficients 
indicating b coefficients of the first independent variable, the second set of P coefficients 
indicating the b coefficient of the second independent variable, etc. 
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Also included in the output data are the number of independent variables N, the 
degree of the model equation P, and the number of input data points (sets) n. Further, 
the significance of the estimated regression equation is indicated by the standard 
deviation, the multiple correlation coefficient, the F ratio, the maximum percent error, 
and the average percent error. 

The values of the dependent variable are also shown for the input data and for 
the computed values obtained from the fitted expression tor each set of data input, as 
well as the residual difference in the input values and the computed values. 

When called for, plotted results can also be obtained which qompare the input 
and computed dependent variables as ordinates to the point interval along the abscissa 
(which is normally one of the physical independent variables, but which can be a unit 
indication of each data set in order of input to program). All computed points are 
connected by straight lines and are plotted with * symbols. The input points are plotted 
with + symbols and the points are not connected with lines. 

Another useful parameter in testing the significance of the regression equation is 
the standard deviation which is estimated by: 


2<Vy,v 

Vi 

S(RES) 

Regression degrees of freedom - 1 


Regression degrees of freedom - 1 

_ 


(35) 


This result stems from the work of Junkin [4] . 

Earlier, the analysis of variance technique was used to test the combined effect of 
the independent variables on the dependent variable using the F statistic. A closely 
related statistic is the multiple correlation coefficient. R. Smille [3] defines this statistic 
as the simple correlation coefficient between the observed values of the dependent 
variable and those estimated by the multiple regression function as given by: 


R = 


S(REG) 

<3 

a Y Y 


(36) 


If the observed and estimated values are completely unrelated. R will be zero and, if they 
are identical the multiple correlation coefficient will be unity. Values in between these 
limits represent different degrees of correlation or the closeness within which the 
regression equation describes the original data. 
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Also of interest is the average of the absolute values of the percent error of the 
dependent variable, considering each fitted observation: 


^ e 1 AVG 


1 

n 


z 



x 100 


During this calculation the maximum absolute error condition can be determined for 
evaluation purposes. 


V. SUMMARY 


A powerful multiple regression technique for Pth degree polynomials with and 
without linear cross products has been developed with which the nonlinear behavior of 
identified independent variables can be related to a given dependent variable. The 
polynomial expression can be of the Pth degree and can incorporate N independent 
variables. The resulting surface fit can be used to summarize trends for a given phenomenon 
and the analytic results provide a mathematical basis for subsequent analysis. 

It is required in the application of this technique that the S matrix be nonsingular 
and that the input data are reasonably well behaved. It is also required that the number 
of data points being fitted be equal to or greater than the number of unknown 
parameters in the model equation. 

Two computer programs have been developed to implement this technique for the 
cases with and without linear cross products. These programs perform matrix operations 
in double precision and evaluate the various unknown constants in the model regression 
equation as well as the standard deviation, the multiple regression coefficient, the F 
statistic, the maximum absolute percent error, and the average of the absolute values of 
the percent error. Further included in these programs is the solution for the identity 
matrix to identify any problems in the original matrix inversion process and a means of 
obtaining machine plots, comparing the computer results to each set of input data. 

Studies of many sets of physical data using the technique without cross products 
have generally yielded excellent results. However, much less experience has been obtained 
using the program with cross products, although test programs have shown very good 
correlation. It is, therefore, felt that the technique developed herein can be highly useful 
to many fields of endeavor. 
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APPENDIX A 


MULTIPLE REGRESSION PROGRAM FOR CASE WITHOUT CROSS PRODUCTS 

A. Input Data 

The regression subroutine is called with the following statement: 

CALL REGRES(Z,YO,LN,N,IP) 

where: 


Z is a two-dimensional array containing the independent variables used. The first 
subscript denotes the independent variable and the second subscript denotes the 
data set. 

YO is an array containing the dependent variables which correspond to the 
independent variables. 

LN is the number of data sets. 

N is the number of independent variables used. 

IP is the degree of the polynomial curve fit. 

These arguments must be assigned a value in a driver program. 

The regression subroutine incorporates several special routines: RDWT, IOWR, 
INVRT, and DMATML. 

RDWT is a general FORTRAN I/O package which is capable of reading and/or 
writing on magnetic tape or drum in FORTRAN or non-FORTRAN format. 

IOWR is an assembly language I/O routine called by RDWT. 

INVRT is a double precision matrix inversion and simultaneous equation solver. 
The call to this routine is: 

CALL INVRT(A,N,M,DETER) 


where: 



A is the input matrix for inversion or augmented matrix for simultaneous 
equations. 


N is the order of the coefficient matrix. 

M=0 for inversion only; otherwise, M is the number of constant vectors. 

DETER is the determinant of the coefficient matrix. 

DMATML is a double precision matrix multiplication routine with transpose 
options. The call to this routine is: 

' CALL DMATML(C,A,B,M,N,K) 

where : 


C is the product of matrix A times matrix B, 

A is the first input matrix, 

B is the second input matrix, 

M is the number of rows of A, 

N is the number of rows of B, and 
K is the number of columns of B. 

Certain variables have dimension changes that vary depending on the number of 
data points (LN), the number of independent variables (N), and the degree of the 
polynomial curve fit (IP). These will be dimensioned as follows: 

DSN(NIP2) , DSHAT(NIP.l), DB(NIP,1), Y(LN), RESID(LN), Z(N,LN), 
SHAT(NIP,NIP), ZBAR(N,IP), YO(LN), BZ(LN), B(N,IP), SVHAT(NIP,1), 
SH(NIP,NIP), XX(LN), PMAX(LN) 


where: 


LN,N,IP are as defined above and NIP is computed as: NIP = N*IP; NIP2 is NIP 
raised to the second Dower. 

SMALLZ, SM ALZ1 , and SMALZ2 should be dimensioned greater than the 
number of data points and the dimension is set in a PARAMETER statement; for 
example: 
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PARAMETER IDIM = 100 


B. Program 

In order to run this program on the UNIVAC 1 108, the deck is set up as tollows. 

@RUN,//T JBNAME,320590,VDARBYBIN406, 01, 100 
@ASG,T 10.F/1/POS/2 
@FOR,IS MAIN, MAIN 

(Source deck for inputting data) 

©FOR, IS REGRES, REGRES 
(Source deck) 

©FOR, IS RDWT, RDWT 
(Source deck) 

@ ASM, IS IOWR, IOWR 
(Source deck) 

@FOR,IS INVRT, INVRT 
(source deck) 

©FOR, IS DM ATM L, DM ATM L 
(Source deck) 

@MAP,I AA, AA 

LIB SYS$*MSFC $. 

©XQT AA 

(Input data) 

©FIN 

as illustrated in Figure 1 . 

A complete program listing for the case without cross products is shown in Table 
2 and a flow chart indicating the operation of this computer program is indicated in 
Figure 2. 


C. Output Data 

The desired coefficients of the model equation are outputted in E notation where 
AO is the computed intercept of the fitted polynomial expression and the B coefficients 
are printed out in ascending order of degree (P) with the first P coefficients indicating 
the b coefficients of the first independent variable, the second set of P coefficients 
indicating the b coefficients of the second independent variable, etc. 
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@FIN 



Figure 1 . Program deck sequence. 
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TABLE 2. PROGRAM LISTING (WITHOUT CROSS PRODUCTS) 
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TABLE 2. (Continued) 
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CL266 • 7 o* ! SE CT = t J- t \ * 1 P-NSEC T * ( K- »> * N *EcT* t 

C.:2«7 7 7* CALL lOHDnK I 2. IUK1T • IbECT ► S* AL* 2 » NtfU$ , I L H r? > 

0:7270 7 st — - GO 313 l»I ,L« • 

02273 79* 300 S u A T I ] R!) .V , I COL > * SH A T ( I WO* , 1 C OLJ *SM A L Z J ( U *S*Al Z 2 < 1> 

C«f27M oO* 1 rfU.’ti ICCL JasHMI iROrt* iCOU 

CC275 cl* 3 1* CONTINUE 

CC277 fa2* IF ( 1 CCL*E«.N.ip> ICOL*0 

CO 30 I uJ* 3*0 UdLalCOLM 

• OCJOM • 6 M * *4 CO CONTINUE ' 

CC30* to* M50 CCHTlfiOE 

Co J iO 06* NKbINOA 

Co 3 1 1 o7 * 12=1 

Q 0 3 I 2 » « • I C L = N * I P 

E D 3 1 3 89* DO 3 3 2 2 i*l»ICL 

GUJ l 6 VC* - DO 1 322 J - I t tflO* ' ' " 

C032 1 Vi* OSf < I 2 » "Shi J , I 1 

C u3 2 2 ' ; M2* i 322 J 2= J 2*1 

C 0 3 2 B I 9 3* H- g 

CUJ2A *■ 9*t* CALL 1 NVKT 1 OSN *NR »Hi DETER } 

Cl- 327 Vb« 00 1 M»Ln 

C U J 3 2 V*** - ' *5t} T8A«0 a Y«AKO*T0 U > ' 

CC33M 9 V* YaA f{ 0 = VbAK02LN 

00335 Vn* - MIP-U*IP*1 

Co 33 6 W* IK0»»1 

CO 3 3 ? 1 00* 00 Mu ]«i|N 

DG3M2 1 u i * !l(l ?li J‘l»lP 

CUJMS , ■ 1-J2* l5ECl-( 1-1 > * lP*NSECl*^U*l >* N SECT*l 

;C3M6 103* CALL I OWD i,R l 2 i l UN I T . t S EC T , S* A L Z Z * N#0 S , \ £ RR ) 

CD 3m 7 1 01* 00 7J J K » 1 * L N 

CO 352 1C>* SVHATUHU»«|l)a5VhATllR07/»n 4 SMALZ2iH)*tT0(KJ-YaAK0l 

CO 3 S 3 le** 7O3 05HATlIROn ( H«SVHAT<lRO**l> 

CO 3b 1> lg/* 710 I MO* B 1 K0«+ 1 

DCOBb 1 0 H » — C - CCHPUTE B’S fOR THE aSSUHEQ HOOEL • • — 

C03&B lUV* C 

CO 3 60 - 1 1 -3 • I R - 1KD7^- 1 

Oj 36J III* CALL .DMATML <DU ►0iU,D3HAT < IN i IR » 11 

0 0 3 6 2 1 ! 2 • 12-3 

CG343 l 1 J* 00 1331 l a |»N 

CO 366 1 I M* — i)0 U31 J°l »lf“ 



u> 

*o 


TABLE 2. (Continued) 


C n 3 7 1 

i 1 b- 

L’ 0 3 7 2 — - ■ 

1 16- 

CG 37b 

t 17* 

CQJ?6 ' 

Uvl* 

Ljsai 

l 1 V • 

00404 

1 2C. • 

CO HO 7 

12 l f 

CCH2 * 

1 2 2 ■ 

CdH 3 

1 2 J* 

C GM 6 

1 24 • 

00 4 2 1 

12b* 

CD424 

1 2 6 4 

C0H24 

12 7* 

C % 4 2 4 

I2n* 

C0H2H 

12** 

C 0 H 3 0 

1 JJ* 

E 0 H 3 3 

13 1* 

C C ** J 4 

1 32 • 

CCJ435 

133* 

D0436 

1 3 "i * 

CO 4 37 

1 Jb* 

"044 l 

1 * 

0C 4 4 l 

i 3 / * 

dJ4 4 1 

l 

0044 1 

l 3 t* 

CC4 1 !? * 

1 4 U • 

L (J 4 4 2 

1 4 1 • 

CC. 44 2 

l 4 2 • 

0044 2 

1 4 J • 

C 0 4 4 2 

I 4 4* 

30 44 3 

1 4b * 

to 4 4 6 

l 4 o ♦ 

0 0 4 4 7 

14/ ■ 

a f J 4 b 0 ; 

, 1 4 3 * 

00'lb2 

1 4V* 

L l. 4 ii 3 

IbO* 

00435 

l b 1 - 

C C 4 b 6 

lb2 « 

Clj4 b 7 

13 3* 

C C 4 6 2 

»t»H* 

00 464 

Ibb* 

0046 b 

l b6 ♦ 

C0466 

1 b 7 * 

JD 4 6 7 

Ibd • 

904 7 0 

Lb?* 

C047b 

1 6 d * 

CDH7S 

L61* 

CDS 74 

J 6 ^ * 

C 9 5 0 0 

J<,3* 

LUbOJ • 

l n *t* 

CO*i 1 2 

1 6b* 

Oc5 l 3 

1 6 a * 

Ovb 1 3 

H/* 

c ;bl 7 

1 6 4 • 

i.jbso 

Uy* 

CCb 2 3 — 

l/J* 


0ll*jl # QBl|2 ( U 

- -- 

AU*TtiAf<0 

DO 23 !°I.N " 

DO 23 J°UI H 

n AO 8 AU'0{I»J^2PAM!>J) * ■ ' 

;; k l T E l 4 ♦ 2 4 > AO 

' 21 TOHMaTI lhU r 3 H A0-»Eie.?l 

DO 7 b y 1 * t . L N 
DO 7bD J»l,« 

OU 7bD *» l » lp 

7 bU D2|I)»BZW>*0CJ|K>**«J*!)**K 

c 

C cnj'.i’UTE THE V *0]hT?i * 

t 

DO 77b 1* 1 ,Ln 
T I 1 J ^ A 0 ♦ 0 * i J > 

YU A lit. 3 YhArtC*T ( I J 

RESIUUI-Y0*|)-TU» 

“ RH»RR*REStOU »«»2 ‘ 

775 CONTINUE 

■ VO AKC a T(UHC / LN 
C 

C COMPUTE THE STANOARO DEVIATION* 

C 

s T AOt v«*5<J1<T l HR/ I LW-N* IP M ' 

C 

C COMMUTE THE lUiL T [ PL t CORRELATION COEFFICIENT ,HHAT i A N D 

c THE F Rat 1 0 » F R A T (0* 

c 

[>0 a^O I a 1 1 Lu 

S4eb=bRt‘C>* ( T j [ j-TOaRC I -■ - 

$ Y * a b V i * t Y o t t I - Y B A H 0 ) • ♦ 1 

dJJ S4ES-SRES+Rty [ D < II • *2 
HHaT = 50ST IS«E0/5YY 1 
tF < KHA T . 6T • J . 3 1 KHAT«2 .U-Kh a T 
F«ATlO*(SH£G/lNMi*-l)»/(SKES/(L«-N*rpjJ 

nitci. * • — 

DO 0 3 LIa2,LN 
*3 X X ( U t I 
L l=M 
L 2*0 

CALL wUlK-3L ( Ll,XL»XfltYBtTT l ‘<g,FL0>wFLDY»-LN>XX, , r} 

CALL ^UlKJL^Z^LiXfltTOfTTiOHiFLDXiFLOriLN^XiYOI ‘ 

WK J T|_ I 6 , lUO^ > STADEVtHHAT.FfiATtO 

I02S FORM A fUHU, I9HSTAN3ARD OE V I A T I ON- E I ** • 9 » 4* » 3 3 H MU L T | PLfc CORRELATION 
■CQEFFJCIENT*E16*?»4X.HhF HATI0*E14,9J 
v/HlTE I 4,22 22 > 

2222 FORMAT t IHJ » COEFFICIENTS') 

■ - - w* I T U 6 « 2 aO" > I < d l ] » J I » J" t ► 1 R > » I- 1 ,H » - ' - ‘ *■*' 

25GD FO^MaTUhO » 5 E2J •*/ I 

arflTEU.uaCl 

.i r. i rc(t>» l u 0 4 J J 

ISOD FCIRMaTU Hy , l6X ,2 ht0i UXiZhTC, 15X ( 5hT0-YO 

DO 2 C w 3 J D I * L N 

... ( ir |T*.i6*2»y2 I TO ( I J » 1 l 1 > »H£S 10 ( 1) ’ 



TABLE 2. (Continued) 


COS30 

CG5J1— 

CD533 

CDS35 
ZCSHO - 
a i ) s si 

C-L'bHJ 
C'jSHM 
C 0 b H 5 
C U S !i 0 
0C‘5SJ 
Ci a b b M 
C 5 5 *■ ■■ 
or:->b7 
Ceb*3 
Cub 6 3 
C-0b69 
OcjS*5 
£0566 


171* 2b2^ kOKMaM iHOt 3e2i*9) 

ZZCL, CCnT t NUE --• - - -- - 

173- i0a^3 r OR*A T < I HU. UX.aHQBSERvEOi l2x»eHC0HPUTE0, l2X t ftMRE5!DUAt I 

- 1 7*M SUPt*C* 

1 7 b • DO Sjl L*i,LN 

- 1 7 o • PMAX»L>"Ai3S l lV0<L)-Y(L))^r(D) 

‘i/7* 531 SUn*IiUM*PHA* IL> 

l‘/d* p£i<A VG = t SOM/LN) * ICC • ~ * 

\ /?• AKAX = PM ax i I ? 

’ U< J * DO 5 J b L » * » L N 

1 a 1 * iF(ANA*-PrtA*<L>> 50 6 i5CS>50S 

1 a ,* ♦ b J 6 A H A A ® H i't A X ( L 1 

163* 535 C0* T J HUE 

- - i 3 1 * • •- • p e i: r, X ” A m A x * 1 c 0 • ' ~ 1 - 

lab* /j * 1 T t ( * i 1 G LS * 1 > PERmAX,PERAvG 

lba* IQ32I FORMaT U H0*22HWAX1m'JM PERcE n 1 E r R 0 R * E 1 6 • 9 i 6 * t 2 2 H A V E R A <* Z PERCENT tFi 

1 6 7* *(?OH b l16»?) 

lau* CALL lUDJOQ 

loV* RETURN 

" 1 9C* E.ND 

L^u OF COMPILATION: NO DIAGNOSTICS. 





TABLE 2. (Continued) 


i*K0«»l5 Ro^r,HOwT 

hsq -oov-oi/i v-d* -( ,Q) - 


subroutine Ki/*T entry point cac 27 o 

keotpk Entry point occ 2 7 3 

KirtiTL« entry point 000320 

— • close — Entry point 0003^2 

uPtlN ENTRY POINT CC 03&1 

iuRo*k Entry point oocjajj 


STyKAOt used; COoEliJ Ow*C‘+02{ D A T A ( G \ 0001271 BLANK C°MM0Nt2I OCOOOO 


EXTERNAL + trtRfcNCEs (BLOCK, NAhE} — 

0003 I Or. R * . ... 

□ <3o*r NEk*2* 

0 C u 6 N I 0 I A 

QZ'J7 ' RI02A •* * • 

C'J 1C NS I OPT 

0311 -- ..... 

SToHrtCit aSSICNMEnI (BLOCK, TYPE, RELATIVE LOCaT|ON, N A hE» 

000 1 00 0 2 R 2 2 iS<5 

0001 ■ 000033 6 L - 
0000 1 000000 1C*0 
0000 -r C00067* iTWt) 

ouoo 1 occuhg kooe 


CD LQO , 1 * 

omoo ’/* 

ColOO 3* 

coioa 

cc lot - s* 
00103 6« 

OfHO** * 7- 

coia^ o* 
co j a* *• 
ClIOA ig* 
LCt ()6 m* 
CC l 10 12- 

LUli 4 * ‘ 13 “ 

C011'* I'i* 

C<j 116 - i B * 

0 LU 2 0 i 6 « 
CC120 Vi* 

C i; I 2 1 lt>* 

Dc 122 - IV* 


C GENERAL PQRTKAN 1/0 PACKAGE E X L C V I L 1 APRIL,IV70 

C HEaO/WP l Tfc M A G TAPE oh drum In FORTRAN OH non-fortran format. 


SUBROUTINE H Dir. T ' - 

DIMENSION aNwayI n»lC^ 0 Cj»*lFL 0 < 29 l.K 0 DEl 6 ».lPACKTlai.KNTwDl 3 J 

DATA I iFLtH 1 ’ » J"l »29>/lHl . 1*2, LH3 ( tH* , |HS. 1**« tH7 f 1*6 . |H? ,2H1Q i 

I 2H t l t 2H 12 i 2H i 3 , 2 h1 I « 2H I& I 2H 1 6 ■ 2HI 7 ,2h) A t2Hl? »2H30i 2*2 1 • 2H22 »2H29 » 

22 H 2 H» 2 K 2 St 2 h ZA, 2 w 27 i 2 N 2 Bi 2 H 2 ?/ - 

0 A T A <K00t.( l ) » 1*1, 6>/O3 0‘+J0CC0 00 CO >0002000000000 i 00 o 1 bogqooooo i 
J 00 0100000 00 00*00 01 lOOOQ 0000 *00 0^0 0 00 CO 0 00/ 

data M/262lHH/al2/lQ737<i|!j24|/ t !3/oO!ODOj/ 

0A7AJBLANK / 6M ■ 


Eh!' <Y KtD I pH 1 1 UNIT »MOOE i l£RR iXtttNIAOS , A«R AY » 

1 OP® MODE £ ■iiKEAO FTN FORMAT* *2 »REaO NO^-FTH, 

t** * * • 

GO Tu 1C 

— • ENTHr AH 1 TEH UUNIT .MODE, 1 EkR,NwOS»A«RAY) 


OCOI Q 0 0 0 7 1 10 L 0000 0 CC 070 lOOF 0001 aOOlt 6 2 aL 0001 OOOIRO 21 L 

Ol'UI CDDJHH 2 2 L cool CGC 2 ll Z**L OCOt 00022 ? 25 L 000 | 00023 ^ 26 L ■ 

OCu 1 OOcajS ? L GCGL 0 C 2 O 52 8 L OQGO 1 000 a 6 l 1 0000 I 000065 IBLanK 

oauu l CObCg 3 1 FL? OCOO CCCItl JNUHS COCO I 00C066 1 qP 0000 [ C0U0H*» IpACKT 

c:.u0 I 000«.o2 I I 00 DO 1 G3CC63 12 OOQO t 00006*+ J3 COOP [ Q0U056 K N T h u 



TABLE 2. (Continued) 


CO 


00124 2u« lOP«*KOUE*2 £ OiaRITe fTu F^RMA** * H , nR I T fc. NO N-F T b « 

Co l 2 5 2|f. K MW U t 1 ) • N ftt> $ * J 1 + I 3 S BUlLO FORTRAN R E C 0 *0 -CO U T rOl « Q R D * 

0 0 12 6* 22* KNTVilD 1 2 ) a J fil ANK H OUMpiy CHECKSUM ERROR ttOHD, 

CO 127— 23* - • <NTU&( 3 >«F.hTwD( n - “ - - . 

0CJ30 2«* GO TO 1C 

Col 30 €*•*•* ' - 

C 0 1 3 1 2** ENTRY CLOSE nUNlT »j0pj 

Lf, 1 33 27 * - tTHO»l ' • 

C £,' 1 3 '♦ 2b* ' GO TO t 6 .6 , 7 , 25 > » JCP fc J0P«H» NO A c T J 0 N , ft E T U R N * 

CCi^ft ' 25* 6 1 T Vj 0 * 2 fi J0K»]t2i V, J T E EOF AND RtwlNQ IUNIT. - 

GDI 3 4* 3v* V [GP a S S "Si AHlTZ E«0-OF-F U.* . 

LC137 31* GO TO 1C 

CC137 32- <.***•" 

LCl^C 23* Cf.TRY OPEM luhlT p JOF J “ “ ■ ■ - 

C C1H2 3H* lFiJOP.GT.il RETURN 6 JOP-2,3* NO ACT 1 ON i RETURN * 

Lul^ • 3S* b [QF*6 k «6» HEMND IUNIT* 

Cc MS 3«* GO 10 1C 

COW?, 3?* <-♦••** 

t(j I H6 3b* ENTRY 1 0*0 V. K U CODE » I UN J T * I A DOR S , ARfi Ay t NfrDS , 1 ERR ) 

C01SC 35* 1M lC0l>E.tfl*ll [CP»H » *H» SEQUENT 1 aL «R[Te. 

Eli; 1 S 2 HO* [ F l 1 tODE.LG* 2 ) IDP°2 « »2» SEQUENTIAL REaD. 

CUISH • Hi* • J P aCK T < 6 » * 1 AcDRS * BEGINNING Oftyn-SECTOK ADDRESS* 

OUISH h2* <.••**« 

OflSS Hi* ■ Jc continue 

CLMS6 HH* l F A C K T < 1 > = t F L D ( I U N t T ► 6 IUN]T IN FIELD-OMA, 

CCI37 "■* V3* l p aCk T 1 2 > = I Bl ANK - 

CL16D S6* 1PaCKT(31»0 

Coibt • **?• 11 IPACNT (H J*kOdE t TOP) til OPERATION code TO IC^1» 

CC162 Htf* GO To 12(3,21 ,20i2 1 »22 i 221 , I0p 

Cn t 6 3 H5* 2 c 1C»<U(1)« 1*1 l*LOC *KNT*Oi 1 > J 0 FIRST CONTROL ACCESS » G H 0 • 

OgH>H ‘DIAGNOSTIC* THE VARtAbLC ARRAY IS OIMENSJOnED* 

CcUR feu* JC*D(2 ) 'NftOS* I1*L0C< ARRAY \ b SECOND CONTROL ACCESS H0«0« ' 

Gl’ 165 SI* Tr*."i<3>" 2 • j 1 *1.0C 1 NNT&D I 2 J J G TH1&D C°nTkOL ftCCtSS POKo* 

Cuibfe » 52* JPaCnM 5|- 3* I 1 *LOC < J CbD l L ) 1 5 NHBR AND LOC OF 1ST C/A POHO* 

GJ167 ’ bJ* GO TO 22 

Go WO -DIAGNOSTIC* THE VaR l A bt^ ARRAY IS DIMENSIONED* 

CCJ7CI S*i* 21 IP A CkT1Sj*N * US*11 + L0C(AHRaY) « NMBR AND LOC OF OaTa *0«US. 

00171 3b* £ 2 CALL U)*YK < ] PACK T I k ROUTINE To ISSUE E^EC REQUEST# — 

v2 1 7 2 5e« IF ( I UP » £ GJ ♦ 6 ) RETURN fi RETURN AFTER RE*IrU 1 NG IUNIT* 

Z C 1 7 R b 7 * ]F<IOP*E^*S j GO TO <2&.&*i[ThO 6 11NQ«|, RETURN, ■2 l «evviND IUNIT» 

001 76 5o* J E RK* 1PACKT 1 H I / 1 2 G ERROR STATUS IN S* pORTlON, 

C y l 7 7 S^* - IF ( t L kN »E u*C) GO TO 2H d NOKRAl O^EKaT |0n INU EkRoR* 

DC 20 1 6 0* IF (lEKH*<3t«^> f»0 TO 26 K TERMINATING ERR0 h * 

00203*.-- 61*- - lERR-2 • C EnU-OF-F lit OR ENo-o F -TAPE f 

CC20H 62* RETURN 

-DC20S - 63* 2 H NwaKnirwOin/jl <5 FROM Ftn CONTROL T(ORD* 

OC 2 06 e>H* IF<HGDC«Eih*I)GDTo2S 

SC? 2 10 6S* na«1PaCKTI*(J-IPACKT1 j I)/II»| 1 6 N* FROh lOPtit 

002 11 «jf»* 2 &ienR*1 

- u02 l 2 bf* REIUkU - 

C 0 2 1 3 orl* 26 WRlTH6.10C J ( 1 PACKT ( I 1 . 1 "1 i 6 I 

00221-— 6 V * ICO FOF- M A T«///5Ch******EaECUT iON TERMINATED IN 'RD»T» \/Q ROUTINE* // - - - 

C* J 2 2 1 /w» U^H 1/0 PACKET IS 6016 1 

002 2 2 •*-: 7 1* STi“<F 

Lli223 72* END 



u> 

ON 


TABLE 2. (Continued) 


U F CUMPiLATfON; 2 OlAGNOSTjcSt 

ASM * AS lOi.H t IUrt« 

asm 1 1 F-o 1 / i v-o i :oo« i-.n > -- - 


1 . 









*11) 

A*R» 


2. 

• Ul 

osooac 

27 

00 


1 3 

0 

ooaooo 

lOftR* 

L 

ao.o,xi i 

3* U 


G 0 0 0 " 1 

72 

1 1 

oc 

OQ 

0 

O^QOCO 


Z* 

io*s 

- -4 . 



" 000002 

74 

04 

00 

13 

0 

COOCOi 



2 • X 1 1 

5 . 


• l u* s 








£N0 



END A^M ERRORS 5 NONE 



TABLE 2. (Continued) 


* V 0 H i 1 S I J* V H T , I N V A T 

HSO-OC*-Ol / 19*0 l 100-4 ,D) — “ 

SUBROUTINE IhVHT EnTHY POUT C00373 

5T.gNA<5L UStlli _COo£.(l l CCOiiii 0 AT A 13 > 0003341 BLANK CQmHONU) 00300(3 


U> 

-J 


EXTERNAL REFERENCES <«L0CK« NAME) 

OQUJ MERR3* 


STORAGE ASSIGNMENT (BLOCK* TYPE. RELATIVE LOCATION, NAME) 


CCU1 

0000 i 6 

11 U 

3001 

33D1H 4 ) 

llSL 

0001 

ooooz 7 

1 1 7 G 

0001 

00 003 * 

123G 

0001 

OOOOH2 

13U 

0001 

03011* 

16*G 

003 1 

03 32 6 0 

1 S 1 L 

3001 

QCC 16* 

2d tG 

oaot 

000206 

2q7G 

0001 

000230 

2 2 0s 

O'Jul - 

- OOQ2A6 

23M G 

000 1 

00033* 

23SL 

C001 

OOC312 

245G 

■ ooai 

000346 

2S0L 

0001 — 

000071 

7SL^ 

53 Jl 

300073 

7oL 

■ ncoi 

QOO 1 t 0 

V 3L 

GOOD 

u ooccu* 

AMAX 

0000 1 

t 000273 

1 

0000 1 

00027* 

U 

DdOO 

1 0032 7 H 

ItiO 

J 0 03 1 

c asaioo 

INDEX 

ccoo 

l 000300 

I N 0 2 

ooao 

0 0 0 3 D 6 

InOPS 

OUOO J 

oaoooM 

IP 1 V 

OuOO 

OOuO 

I 0 0 J 2 7 S 

0 JOOUJO 

In 

s i on 

0 0 30 

I 030273 

J 

□ 020 

l 0C327* 

K 

0000 

[ 000277 

L 

0000 I 

00027 1 

nn 


00)01 )♦ SUBROUTINE 1 N V R T t A , N , M , D E T E ft ) 

00103 2* PARAMETER |U)M a 60 

50103 3* C (1ATR]A INVERSION AND SIMULTANEOUS EQUATIONS SOLYER 

GolU3"~~4* L A^JNHUT WAT*!* FUN lNVE«lSON OR AUGMENTED MATrU FOR 5 JME* ERS* - -1CS0003C- 

grifi? s* C OP COEFFICIENT NATRU I^SOdORQ 

C0I03 ' - 6* L M*3 FOR tNVtrtJsON ONLY 

001 03 7* C M**nUKBER OF C0N5TANT VECTORS 

c 0 l 0 3 • • U* C deter* determinant O f coefficient matrix 

00104 V * DOUrtue PRECISION A ( t f , DC T E K , S I GN , AM A X 

Ol/LflS- 10* DIMENSION l^iVllDlM>,INOEXlIOlM,2> - 

GOUJ* 11* DETKrfni*CijQ 

CL‘107 - 12* SltiN-J.OOU '* ‘ 

Cat i a 13* d a 2 o j * 1 1 n 

ooi i3 - in* • - 2a rpivi jj«o 
COILS IS- NN»N*li 

Cal 16- 1 ft 4 — ■■ DO Id* K"l |N ■ - - - 

Go i 2 1 ) 1 • AMA X«U-ODO 

00 1 22 - la* MJ 00 76 I ■ I . N 

C0125 IV* \f 1^^1)1-1)50,76^0 

CO i 30 2 0* 53 DU ?5 J 4 liN 

0C133 21* IF l IPIV{J)-1)55,7S»2S0 

CC136 2 ?.+ - SS I MO" l J- 1 ) *N* I - “■“* 

C0J37 23* iF<ArtAA-DABSlA<lNO))) 6Q»7S>7S 

COM2 - • 21* 6s !!»' I 

Cl) L H 3 2S* IC*J. 

GOIHH -- 2o* *MA AdSl A < IN0) > 

CO)MS 27* 73 CONTINUE 

CC 1 *7 — J a • 7 6 CONTINUE 



OJ 

00 


TABLE 2. (Continued) 


CaiSl 2 9 * lP]Vt]C»-lPlVflC»+l 

CD 1 b-2 3- 4 .. [F l IK-tCl 93»nS > ?5- 

CD1S5 • »] AGNOST Jc* CONTROL CAN NEVER REACH The neat statement 


C 0 1 5 5 - 

J J * 


IF ( IN- 1C \ ?Q , J 15 »90 

CO 160 

32* 

90 

S I « N B - S IGH 

cum 

33* 


no UO L« i ,Nn 

Co l 6 a 

Si* 


I HU* l L- t ) *N*J R 

CC16S — 

3 ‘.j * 


I Nu2» ( L - 1 1 *N* j c 

0 1; l 6 6 

3o* 


A H A X » A { [ NU ) 

£0167 

3 7 * 


AI1«U) = M1NU2I 

00 170 

3 l » * 

I 1 J 

A l I N 0 2 1 * A .1 A X 

CO l 7 2 

37* 

1 lb 

I \'|)E A ( K » 1 ) * I r< 

Or, 1 7 3 

h.;» 


jh.;EX (K,2)«tc 

t ’M 7* 

•1 i * 


IN.,* 1 IC - f i *N+ jc 

C f j 1 7 S 

42* 


A f ! A A 3 A t ! ND ) 

CO l 7 4 

4 J • 


l>£ T £ *= 0 r. 1 L!(* AH AX 

CO 1 7 7 

4 *j * 

140 

A(]NUJ*t.oDQ 

rc^uo 

4 b * 


0 0 l bC Lai ,H« 

C 0203 

Ao • 


1 jJ i) * l L- l ) • N * J C 

co 2 n*» — 

4 7 * 

ISj 

A l I N u » 3 A f 1 r< 0 i / A M A A 

13 & 2 0 6 

H < * 


no 1 '3 i L * 1 1 N 

Cl' 2 l i 

4 y * 


JF IL-1C fUS, j8l,|6b 

C*J2 1 4 

S j • 

16b 

HIO* (IC-I M fg * L 

Ui)2 l Si 

3 1 • 


A M A * - A < I NU ) 

00 2 1 6 

S2* 


A l IUU l a C ,003 

C02 1 7 

b J * 


do i ao i • i * nn 

CO 2 2 2 

t» A * 


I NO- < I - 1 1 * N 4 L 

C.02Z J 

b!>* 


I N o 2 B l I - I ) * N ♦ I c 

C ;; 2 2 4 

b 6 • 


A i I Nl> )=AUN0)-MJND2) 

Cf7225 

b7* 

1 a j 

CO.Nl J HUE 

: c 2 2 ; 

'jb* 

1 0 1 

continue 

l 2 3 \ • 

by * - • 

1 *2 

co‘;M nue 

0 C 2 3 3 

t>C * 


DO 2 3b m.N 

002 36 

6 l * 


L * N 4 l " I 

c r: 2 3 7 

O 2 * 


1 R* I fiOLX ( L 1 1 ) 

C L- 2 H o 

6 3 * 


1C*1NUEX tL»2) 

C r 2 M 

6 4* 


IF <1R-IO210i23S}2J0 

0:12*14 • 

ob * - 

2 1(J 

DO 2 3 C K» 1 , N 

£,0247 

o ei • 


I NO* ( 1 ff- 1 *•«♦!£ 

CC2S0 

6 y * • 


IN02« IIC-1 J * N * K 

0L2bl 

o ti * 


A M A X = A 1 1 N D > 

C02b? 

6-7* 


A U ND J *A < lhV2 > 

CC253 

73* 

230 

A 1 L Nl) 2 J * A H A X 

C0 21iS — 

7 1 * - 

23b 

continue 

CO 2 b 7 

7 2* 


detlk=si un*oeteh 

00260 

7 3* 


ret UNn 

C 0 2 6 I 

/4* 

2b U 

M«- 1 

Cp2 6 2 

7 j* 

231} 

RE TOkN 

C 0 2 6 3 

7 6* 


end 


i N o OF COMPILATION: 


1 DIAGNOSTICS. 



TABLE 2. (Continued) 


WfQR.lS DrtAFrtL t LiHATm. 

Hgo-of -Ci / i»-oi :ov-i ,o ) 

SUBKuUTlNE OHATMU entrt POINT 0001*0 

._„szoRA£e ..ust 1 1 >... 00017 * j- data 10 ) oooa^ai -blank cQnMoHU>„aoQOoa 

external RertRe Nets ibluck, name) 

03 U 3 NEHHJS 


— S T u rt a G L Ai LAMENT )BLQCK, T T P E , RELATIVE LOCaT[ON, NaMO 


.. COO) 

OCOQiJ U 

0001 

acoo*i 

l 30G 

000 1 

coco* 7 

J 3 5 G 

□ 00 1 

0 DO 1 D 2 

1 

-Oaoj - ODcgiS 

-2L — 

OJUl 

□C00H6 3t 

OOJ] 

003052 

HL 

00 0 0 

D OOOCOO 

CO 

OODO 1 

I goooos 

] A t 

oaoo i ocooo* 

1 A 2 

OOUU 1 

000002 ' 1 A3 

- - OOOC ] 

i C32C0 7 

101 

' cooo 

I oojclc 

I B 2 

0000 ; 

I 0000 1 3 

IB* 

— OCOO-I-OODOO^ 

IM~ 

OGOO 1 

OOOGJ** In 

QCGO 

030022 

INJP* 

coco 

1 ooaois 

LA 

0000 : 

I 000016 

LB 

□000 I 000012 

LC 

• OCUO J 

ooaotl LM - 

OODO ] 

I 0 G 03 I 7 

IN 

0000 

I OOQOlR 

LP 



- 




00)01 

1 ♦ 


00)01 

2* 

C* 

00 10 1 

J* 

c - 

oc 101 

H* 

C 

00 1 0 1 

•j* 

C 

Co 10 J 

— 6 * - 

C 

Co 1 □ 1 

/•- 

C 

00 i U J i 

« * 

L 

G01G) ‘ 

?• 

L 

COLO) — 

) O* 

C 

CJlCN 

1 1 * 

L 

00 1 0 1 

- 12*-“ 

c - 

Cl 0 1 0 1 

1 J* 

c 

aa to 1 — 

M* 

C 

CO l 0 ] 

lb* 

c 

OclOl 

)6* 

c 

ooioi 

1 /♦ 

c* 

CO 1 1)3““ 

• 1 Cl * " ■ 


(<U 1 1)3 

1 '} * 

c 

CJI03 

20* 

C" 

00 1 0 M 

2 1 * 


OCiUH 

22* 

L 

CO l(]H 

23* 

c* ■ 

00105 

2H« 


00 106 

25* 


ecu? • *■ 

26^ 


LOUC 

4?* 


CO t 1 1 

20* 


00)13 

2V* 


conn 

33* • - 




* i TH TRANSPOSE O p T IONS 


SUBROUTINE 3f1ATMLlCtA»BiM,N 
ABSTRACT 

SEVERAL MaT«IX MULTIPLICATION ROUTINE 
*h£RE, M IS THE NUMBER OF ROnS OF U) 

N|S ThE NUM 3 EH OF *Qft 5 OF <S) 

k i S the number op columns op c a 1 oR ib)t 

T rt f N S P 0 $ £ OPTIONS A C^TROll£D BY ThE 5 I fiNS OF M AND N, 
THE fOLLomlNG PRODUCTS MAy (Jt OBTAJNtO 
M A^D M POSITIVE 
M N£GAT l V E FOR l A ) T 
hi NEGATIVE FUR (BIT 
M AnO N NEGATIVE 

transpose 

H IS ThE NUMBER OF 

is The number of 


C C ) * ( A 1 < e I 
t c > - ( a m a J 
1C j * ( A J u 1 T 
(C ) = l A > T l 9 ) T 
*hE*t T INDICATES 
IF rt IS NEGATIVE, 
IF n [5 NEGATIVE* 


rO*s of 

R OiK 5 Df 


I A IT 
(BIT 


OUTPUT ARGUMENT * C 
dimension C 1 J > 


INPUT ARGUMENTS * A»t)iM*N,K 
DIMENSION A>1)|BI|) 


double precision co>c»A»d 
] A3* 1 

IM* UdS IM J 
I M» l A 3 S (M J 
I F I M • LT ♦ O f GO TO 1 
I A 1 a l H 

t A 2 " l ■ • ■ ' • 


U> 

VO 



o 


TABLE 2. (Continued) 


0 011 3 3 1 * 

-C »j 1 1 6 3 2* l 

00117 33* 

00120 ■■ -2 
C 0 i 2 2 3 b* 

0 0 12 3 3 ft * 

C 0 1 2 H 3 7* 

0 p 1 2 S * 3d* J 

DC 126 3 Vi* 

CC127 Hj* i 

0 0 13 2 MM 

C i; 1 3 i M2* 

C0 13M *4 3* 

00 137 - MM * 

Gill MO MS* 

CO J H \ Mft* 

C 0 1 H l M 7 * 

C 0 M s M3* 

COJM6 MV* 

0O1M7 SO* S 

CO 1 * 1 S 1 • 

0 G|02 b2* 

OUISJ S3* 6 

CIO I S S " ’ bi* 7 

0 u 1 S 7 S ** 

00 16 0 5 c* • 


00 TO 2 

-- I A l • a. \ 

IA2MN 

1 f 1 N ‘ « tT # 0 1 GO TO 3 

1 0 1 55 I 
182" IN 

GO TO M 

t'JlH ' - 

M2* I 

DO 7 LM* | i 1 M 
IC*LH 
1 33M 

DO 6 L P * 1 i K 
CO ■ 0-200 ' 

L A a 1 A 3 
L rt a .1 b 3 

DU * LM a I ♦ IN 

CO ■ CO + A<LA|*tf<tB) 

LA^LA*JAi 

L H a L U ♦ I 3 1 - ■ - ■ 

C < lC > * cu 

tc * LC * JM 
I6 3 s ld3*Ii32 r 
I A3«M A3* I A2 ' 

RETURN 
E ?J 0 


e^u ok compilation 


NO 


diagnostic^ 



TABLE 2. (Continued) 


■MiAPilS Aa>aA 

*01 

l* Lie syss*mSfcs» 


AuORC^S, L4n I TS ... OiilGa; 033207. Q400QC 054642 - 

STAhTIM* AUU^ESa 033 167 

AUKOS UECIMaL 1344a IbANK 6563 (J3ANK 


SEGMENT main .001000 0 33 2D 7 


NbrtTC*/rOH 

i 

00 J 0 0 0 

501321 


NkHL«*'MSFC35 

1 

00*022 

051*10 

0 

«t APA*/H’3r tS 7 

* 1 

GO* m 

051 3o3 

2 

ALOOS/Fo'^bl 

1 

00 1 304 

03*421 

2 

CS I G6v/ >C4 u20 


00 1 H 22 

0C1 7 *0 

0 

cewMtfn/bC4;:2Q 

1 

GO* 7 11 

00 1740 

0 

2 






0LA3LV ^ St 4 020 


00*741 

U JJ 

U 

y 

fr Ok^O 


00*004 

003063 

c 

1 

N*E FS/MSfCbS 

i 

30*064 

CO 3 3 2 0 

2 

fcf TCMi/FOi<b? 


00^321 

CC362Q 

2 

CYK00V/SC4C20 

1 

003621 

353*46 

a 

2 

0 

CaCCH Y/5C4070 

L 

JG 3 647 

003*70 

CXKOO V/SC4Q20 

1 

UG 3 67 1 

P03 7 l 6 

C 

Q 

cokcaf/mSfc 

l 

DO 37 l 7 

0C4C 70 

a 

StT INT/3C4020 

1— 

u0 H U7 t 

0&412J 

0 

ChOLL V/SC 4020 

1 

004 ] 2 * 

004206 

0 

5 

Ct(ONLtt/3C4 320 

1 

00^207 

0046 1 4 

c 

D 

* 

CLIMSV^SC4020 

— l ... 

OC^MS 

C0S470 

£ 

0 


3 

GG<1 


2 

C YSCL Y/504020 


00*471 

CCS*73 

0 


. 3 

GG<j 


2 

CXSCLV/SC4C20 

1 

03*6/4 

QC&C76 

0 


3 

GG* 


2 

tt«NLv/SC4020 

- ' 1 

00*077 

00*263 

0 


DH0OQO 054*42 


040000 04000! 
040002 04DOS3 
0 4 C U 5 4 0 4 U l 14 
040115 C40lb3 

BLANk^COMfiON 
C 4 Q 1 5 4 040*67 

bla^k^cummOn 
CH017C 0^^276 
BL AM^sCOhMQN 
04027? 040310 
0403 1 l 0^03 3 1 
OH0332 040367 
0 4 0 3 7 0 li 4 0 3 7 6 
ttLAfvKSCQtlHOr^ 
0403 77 0 4 0 4 fl 7 
a L A N K $ C 0 M M 0 N 
G 4 0 4 10 0 4 Q 4 16 
BL AnKsCOHMOn 
Q40417 040^40 
040441 G4Q446 
BlAkksCOMmOivI 
040447 C40 1 *^3 
BLAfcKSCOwMOnl 
0411464 04Q s 3b 
BL ANAJCOMkO* 
0405J6 0 4q65 I 
BLAliHsCOrtMOw 

C406B2 O40 ? SO 

8LAKK*C0«h0u 
040/Ul 040727 
BL ANrfsCQUfcOH 
C40730 C4c 7 42 



TABLE 2 


CEaRLSi/-SC*320 

CSETCWSC402Q 

Ci£TMV/SC402Q 

OAXtS^CHCZO 

VCH AF?v/bCH020 

R 1 TE2 ;/SC4020 
BHL0TK/5C4G2Q 

o PMP/SCH02C — 

CC AhK A/SC*iu2u 

Tj»BU V/SCH22Q 
NdDC Vl/F 0*67 

*► TV*/F0* ' 

*CLCS*/i4iM-C57 
M,\PLK» / m3kC57 
Nt»$dL*/F OR 
NuPOAlfr/fOR 
NtifCOa»^F°rt 

T 0*57 

K [ N I N :* f M S r C 5 5 
CPLQTV/$C4d2D 

CLINEv/SC*4020 

Vi»CtVl/i»C«*«<l — 


X$Ct-V)/5t'»J 20 
C aPIO T/ bC*C23 

C 4Prt^ */5C4 223 

CpK*Tv/SC4!?20 • 

C&RD J V/5CMO20 
COAOY+S SCH020 - 

Cti4iT*/SC402G 
C*A*G:g/iiC402a 

CmjLNS/SC^vZO 

bKoy/rtif c 

Trace - 

ir-Cf RH/HSKC 


2 

.--1 00*25* 036363 ... o .... 


2 


• 1 

00*364 

DC6 4 2 H 

0 




2 

’ 1 - 

00*425 

006503 

0 




2 

... J.... 

- 00*504 

006 7 07 

0 




2 

1 - 

0 D 6 7 10 

307153 

C 

j 

0 0 7 154 

007404 

0 

1 

00 7 435 

037466 

0 




2 

- ! 

0C 7 4*6 

007/35 

0 

3 

GG<J 


2 

1 

00 7 706 

007/57 

a 

3 

G&G 


2 




c 

1 

0 D 7 7 6 0 

310M3 

2 

l 

ClG 1 1 4 

3 1 0 l 36 ’ - 

— 

t 

01-13/ 

C 10305 

2 

t 

010306 

□ 104/4 

0 

1 

0 10475 

013532 


1 

0 J 0533 

013&66 





2 

| 

0 1 0567 

011321 

2 

1 

0 1 1022 

01 1*63 

2 

1 

0 t 1254 

01*436 

0 

J 

\jijG 


2 

1 

011437 

3)2:22 

U 




2 

- 1 

0 l 2 a 2 3 

01 2 141 

a 

3 

GG 5 


2 

1 

U 1 2 | 4Z 

012257 

a 

3 

OG<j 


2 

1 

a 1^260 

0 1 2 5 3 3 

(i 

3 

ggg 


2 

1 - 

0 1 2534 

t>l 2* 1 * 

0 

1 

0 126 17 

0132 10 

2 

0 

3 

GGG 


2 

I 

U 1 -3 2 1 1 

014 l l 5 

0 

3 

GG* 


2 

• - J • 

•014116 

014 7 l 0 

0 




2 

1 

014711 

015312 

0 

3 

ggg 


2 

l 

0 l 5 0 1 3 

015071 

0 




2 

_ i 

Ci 1 507 2 

015141 

CJ 




2 

1 

0 t * I 4 2 

015222 

a 




2 

1 

il 1 5 2 2 3 

0 16-346 

u 


2 

0 


. (Continued) 


BL ANkf COMiiOM 
0 HC 7 4 3 04Q760 - 

blankscokmQn 

340/61 0*G77l 

tJL*fiK«C.OMhON ' 
040772 G41G06 
6L ANksCOHmOn 

041QO7 041042-- 
•iU*NK«COMH0N 
OH 1 0 4 3 C *t 1 0 6 2 
QH 1 06 3 04 1 I 1 l 
OH l J l 2 Q41H7 1 
BlAnk *C0h«0n 
341472 CM]S5o 
BL AnkSCOMhOm 
□H iS57 CHI S/£, 
blank sCommOn 
04 1 b 7 1 0 h 2 l IJ 
0 H 2 1 L 2 042 Tb! 


342 (52 042177 
042200 042201 


04*202 0 H 4 HQ 3 
04 H Hij4 044473 
044M74 Q 44 5 I S 
044616 0445*4 
OlANKSCOMhOft 

04HbS5 <144*33 
a L A N * * C 0 fl h 0 N 
04HA34 044*44 
BLANKsCOMnON 
044445 044*55 
BLANAStOtlMON 
04 H 6 b 6 0 4 4 7 J | 
BL ANKiCOMflO^ 

C 4 4 7 3 2 04476?) 

8LAj,KSC04hOn 

044751 045005 
SLANKsCOmmOn 
OH^OO* 046 l 2 J 

BLAfiKSCUrtflOrj 
04&122 C4520i 
BLANK SC 0>1MO|* 
04*207 045222 
fi L ft NK SCll rip ON 
U4 5223 046264 
BLAnk sCUHhOn 
CH52o5 045300 

uLAjNKjCOHhON 

C4530 1‘ 0 4&3 1 j 
8 l ank * common 

OH5j 12 u 4 6 3 1 7 
UH*320 046402 
04^403 i) *t 5 4 5 6 



TABLE 2, (Concluded) 


CTOaCu/llSKC 

1 

01&3H7 

81S*4S 

0 

045457 

045*30 

aoXLT*/CSC- 

) 

•015646 

0 1 6500- - 

- 0 

- 045631 

045*47 

Cli>£UT/SC4C2a 

1 

0 1 *5C 1 

C J 7*53 

0 

0 4 5 0 5 0 

046033 

* - * - 

3 

CG r * 


2 

OL ANKsCUM mOn 


1 

U l 7 6 5 4 

S2C 2?4 

2 

04*0 34 

C 46044 

ti 0 U T s / 11 5 f L 5 7 

l 

U202SS 

02 1 i 66 

2 

0 4*045 

C46076 

NFtfTWtf:>FCb7 

1 

3 2 1 16 7 

U22C7H 

2 

C 4 * 0 7 7 

C46 1 1 5 

N I OC*>2/ mSf C 5 7 — 


02 20 7 5 

022^54 

z 

C4* 11 6 

046240 

NFCHKj>/KSfcB7 

1 

022255 

023073 

2 

04*243 

C 4 6 H 15 





4 

C 4 * 4 l 6 

ti 4 6 4 * 7 

NTA6S/M$FCb5 




2 

0 4 6 4 7 c 

046556 

N l INPS/HSFCB7 

1 

02 3 U74 

0 2 4 * 4 0 

2 

04*557 

0467 7C 

ByHTi/F UftbB 

l 

02^64 1 

C Z 4 7 r c 

2 

04*7/1 

0 4 7 C U 2 

NeXP55/F0R37 


32‘ 1 7C 1 

C 2 '» 7 6 4 

2 

C 4 7 U U 3 

0 4 7 0 1 2 

G^G <CUhMOh UlOC*) 





Q4?ei 3 

0 4 7 1 4 6 

C{<U]KL/5CHC?C 

1 

02*765 

025*02 

C 

0 4 7 14 7 

04723 t 


J 



l 

blankscgmmOn 

IUFNT/5CHC2G 

1 

02^33 

024*45 

c 

0* 7 Z 32 

051027 

HiMONlTQR/nBFCSb 

1 

02*640 

0277C.S 

2 

051330 

0 5 1 6 J H 

NIERS/F0HS2 

j 

02 7 75A 

C3CC4C 

2 

051615 

051744 

NUtilJF"S/FOi<‘i l 


0 3 C C « l 

C 3 C I GO 




tHUS/ilbF CS5 







NEKSS/FORb / 

1 

030 10 : 

03cM5 

2 

051745 

052121 

b L A N K S C 0 M M 0 N ( C U M M 0 N 

BLOCK J 






D H A T H L 

1 

03C*436 

D3c6 J3 

0 

052 1 32 

052174 

. . . .... _ • - 

— * - 



2 

13L *NK SCOMrtON 

IrwRT 

1 

0 3C o 3 4 

0 3 1 2 S 1 

0 

052 1 75 

052530 






2 

bLANKSCOMHON 

I 0 ‘A R 

1 

. 0 3 1 25 2 

03 i2$4 




hU Vi T 

1 

03 • 255 

C3 1 *56 

0 

052531 

052*57 





2 

BLANKSCOMhON 

RtGRE* 


03 1 45 7 

033 1 66 

- 0 

0 5 2 0 6 C 

054^41 





2 

BlANK%COM^ON 

MAIN' 


333167 

C 3 3 2 c 7 

0 

05H542 

C 5 4 6 4 2 


2 BLANKSCOMmON 


„ S Y ss *kL| 8s ♦ .LEtftlL. M 5 7 ~ 5 . . 

END OF' COLLECT I UN - T I M L 3*518 SECONDS 

ii X li T A A 
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DIMENSION Z( 3^ 10),Y0C 10) 



DIMENSIONED 


SYMBOL 

STORAGES 

SYMBOL 

STORAGES 

SYMBOL 

Y 

10 

fi£$ID 

10 

Z 

ZBAR 

3,3 

YO 

10 

BZ 

SMALZ1 

10 IM 

SMALZ2 

id in 

SH 

PMAX 

10 

FLDX 

12 

FLDY 



VARIABLES 


STORAGES 

SYMBOL 

STORAGES 

SYMBOL 

STORAGES 

3,10 

SMALLZ 

lorn 

SHAT 

9,9 

10 

B 

3.3 

SVHAT 

9,1 

9,9 

XX 

10 

ADARY 

22 


12 


Figure 2. Operation of the computer program (noncross product). 






SUBROUTINE REGRESf Z,YD # LN,N,1P ) 
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* 



-P* 

KJ\ 


Figure 2. (Continued) 
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SUBROUTINE REGRES( Z,VO,tN,N,TP) 




& 


K 


TRANSFER TO SUBROUTINE 
IORDWR 

lUNIT, I SECT, 5NALZZ,NUOS, TERR 


\ \ REPEAT TO IOO \ - 


SVHATt !R0N,1 )=SVHAT( 1R0UI, I HSMALZa KM V0< K>-VBARO> 


Figure 2 . (Continued) 










SUBROUTINE REGRES< Z,VO # LN,N, IP ) 



4 ^ 

<1 


Figure 2. (Continued) 


























































GENERAL FORTRAN I/O PACKAGE EXEC VIII APRIL, 1970 

READ/WRITE PIAG TAPE OR DRUM IN FORTRAN OR NON-FORTRAN FORMAT. 


DIMENSIONED VARIABLES 


SYMBOL 

STORAGES 

SYMBOL 

STORAGES 

SYMBOL 

STORAGES 

SYMBOL 

STORAGES 

SVMBOL 

STORAGES 

ARRAY 

1 

ICWO 

3 

IFLD 

29 

KODE 

6 

IPACKT 

8 

NTUD 

3 










VO 


Figure 2. (Continued) 



Vl SUBROUTINE POUT 

o 




DATA! IFLD( I >, I = 1 , 29 >/ 1 H 1 „ 1 H 2 , 1 H 3 , 1 H<*,IH 5 , IH 6 , 1 H 7 , 1 H 8 , 1 H 9 , 2 H 10 , 2 H It, 2 HJ 2 , 2 H 13, 2 
H 1 H, 2 H 15 , 2 Hl 6 , 2 Hl 7 , 2 H 18 , 2 H J 9 , 2 H 20 , 2 H 21 , 2 H 22 , 2 H 23 , 2 H 24 , 2 H 25 , 2 H 26 , 2 H 27 J 2 H 26 1 2 H 29 


[ START 




DATA! KOOEt I ) , I - 1 , 6 1 /OOO^f 300000000 , 0002000000000 , 000 1500000000 , 0001 000000000, 00 
01100000000 , 000 ^ 000000000 / 












H 


GO TO 10 


!TW0=I 


cm 


GDT0< 6,6,r ( 2 S ),JDPcO a p,‘,,NDACTrON,RE™ m . - HlTWrWOP-l.Z.MHlTCErFAHPHEMimiUWIT. 


hhl 


* IOP=5«=5,WRITEEN9-OF-FJLE. 


GO TO 10 


JE. 


JOP.GT* I 


VES. 


RETURN«J0P=2, 3,NDACTI 0N,RETURN . 


X 


. m 


j£ 


|QP=t«=t,fiEMIWD[UWIT. | ) GO TO 10 ICOOE.EQ.l JEij JOP^^SEQUENTIALUIRITC . 


lL 


JE 


J ICO0E.EQ.2 


I0P=2«-2, SEGUENTIALREAD , 


M4 


IPACKTC 6) = rADDR5«0EGlNNINGDRUn-*SECTt)RfiDDRE5S. 


f 10 I 


CONTINUE 

IPACKT< 1)=IFLD( JUNIT)«rUNlTlNFIELD-DATA . 
IPACKT(2)=IBLANK 
IPACKT(3)=0 


11 


1 PAGKT( M ) « KOO£< i O P O«OP£ftATl(WCOO£W€< 


3 » 




Figure 2. (Continued) 
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D 1 B E HI S I 0 N E 0 VARIABLES 

SYMBOL STORAGES SYMBOL STORAGES SYMBOL STORAGES SYMBOL STORAGES SYMBOL 

IPIV IDIM INDEX IDIM.2 


Figure 2. (Continued) 


STORAGES 


7 



SUBROUTINE JNVRT< A,N,fl,DE7ER) 




Os) 


Figure 2. (Continued) 
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SUBROUTINE INVRT< A,N,n,DETEf?) 



Figure 2. (Continued) 
























DIDENSI ONEO 


VARIABLES 


SYMBOL 

STORAGES 

SYMBOL 

STORAGES 

SYMBOL 

STORAGES 
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Figure 2. (Concluded) 















Also included in the output are the number of independent variables N, the 
degree of the model equation P, and the number of input data points (sets) n denoted in 
the output as LN. Further, the significance of the estimated regression equation is 
indicated by the standard deviation, the multiple correlation coefficient, the F ratio, the 
maximum percent error, and the average percent error. 

The values of the dependent variable are also shown for the input data and for 
the computed values obtained from the fitted expression for each set of data input, as 
well as for the residual difference in the input values and the computed values. 

Plotted results are also obtained which compare the input and the computed 
dependent variables as ordinates to the point interval along the abscissa (which is 
normally one of the physical independent variables, but which can be a unit indication of 
each data set in order of input to the program). The computed points are connected by 
straight lines and are plotted with * symbols. The input points are plotted with + 
symbols and these points are not connected with lines. 


D. Illustrative Problem 

For the purpose of illustrating the use of the program and its capabilities the 
following arbitrary equation was selected for use: 


Y = 9 - 6Xj + 7Xj 2 + 5Xj 3 

- 4X ? + 1X ? 2 + 7X 2 3 (37) 

- 6X 3 + 7X 3 2 - 6X 3 3 


From equation (37) the following set of input data (Table 3) was developed 
containing exact dependent variable values (to three decimal places) for arbitrary values of 
the three independent variables. 
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TABLE 3. INPUT DATA, DEPENDENT VARIABLE VALUES 


Data Point Number 

Y 

X 1 

X 2 

X 3 

i 

5.719 

mm 



2 

7.111 


0.3 


3 

5.216 

0.2 

0.6 

0.6 

4 

5.592 

0.2 

0.6 

0.5 

5 

10.915 

0.9 

0.0 

0.5 

6 

7.569 

0.7 

0.0 

0.6 

7 

7.450 

0.3 

0.9 

0.7 

8 

9.618 ] 

0.9 

0.3 

0.6 

9 

5.019 

0.5 

0.6 

0.7 

10 

8.021 

0.3 

0.8 

0.2 


The input for this sample problem was read in through namelist NAM 1 as 
follows; 

Col. 2 

SNAM1 

Z = 0.5, 0.4, ..., 0.2, 

YO = 5.719, 7.1 11, ..., 8.021, 

LN f 10, 

N = 3, 

IP = 3, 

S 

The computed results are shown in Table 4 for this example problem and the 
plotted results are shown in Figure 3. 

The plot graphs the point intervals on the X axis against the exact and computed 
dependent variables. The exact dependent variables are plotted with (+) and the points 
are not connected with lines. The computed dependent variables are plotted with (*) and 
the points are connected by straight line segments. 

Inspection of the digital and plotted results reveals that the computed regression 
relation yields an almost exact representation of the input data for this arbitrary case. 
Experience with many sets of physical data have also shown excellent results. 
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TABLE 4. OUTPUT DATA LISTING FOR ILLUSTRATIVE PROBLEM (CASE WITHOUT CROSS PRODUCTS) 


J jis J L f4a 1C 

a U * •09/989665+21 

SlAKU«KD ot» I A T 10'ia . I H 9 3 7 3 S6 8 -0 H MULTIPLE CCHStL*TlUN C OtP F 1 C I E NT » . 9 9 V 9 99 97U*00 


0 COtfF 1C I ENTS 


-•bV-' + bzc/z+ai 

.699366400*01 

• 7 jaaOVHb9 + J 1 

-•599687073^01 

Tu 

YC 


computed 

• r* 7 1 6 1 V 9 ? H + 2 X 

• 5 ; 1900227*01 

•711 lC^UOO+01 

i 

• 7l l 1CC20C + CI 

m\, 2 + m 1 

. 52 1 &CC l C3 *C 1 

• 3 :i 9 2 j w 0 a 1 + G 1 

« 5^9 199602*0 1 


• iov i5r.oe« + c2 

• /bO 9*000 j + i | 

• 756B9V226*G 1 

. I ♦, 4 V V 7 Q ? V ♦ C I 

. 745CCC2 A'?*C I 

• ,,, 6i6C £ 0 C 3 + 2 1 

• 9 h 1 7V935 f? *Dl 

♦ ti Cl l h 9 •/ 9 9 9 ♦ y J 

•5C19CC327+E1 

• ’jl21u'J 003 + 01 

• GC2099l>7h*0 1 


.5000*0704+0 I - • SOOCO 1 669 + G 1 

• A996 0 H674+Q \ -• 599740306+0 1 

vo-YC , * - 

RES I DUAL 
-.2304 1 6*79-05 

-.2eo 

-• 101327696-05 

. 1 19209^90-05 

-*SH63ti59&6-C5 - 

•77^660362-05 
2205^7 1B6-C5 
. * 6 H 3 7 3 0 i6^-0S 

--327o2bb46-C5 

• 429 153 * 1 ^ 2-05 


F RAT 10“ ♦ 1 904S029 1 ♦ i I 

.99*636566+00 


i-jfcZ Jr or- r» t, - Ct. i< 


L *< U P B 


1 0237299C-C3 


AVERAGE percent LKHOr 


&17U97b60-Ofl . 



NOTE: THE INPUT DATA {+> AND THE COMPUTED RESULTS (*) 

FALL AT THE SAME POINT FOR EACH CASE. 

I I I l I I I I I I I 

0 1 23456789 10 11 

POINT INTERVALS 

Figure 3. Comparison of input and computed values for illustrative 
problem without cross products. 
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APPENDIX B 


MULTIPLE REGRESSION PROGRAM FOR CASE WITH CROSS PRODUCTS 

A. Input Data 

The regression subroutine is called with the following statement: 

CALL REGRES(Z, YO, LN, N, IP) 


where : 


Z is a two-dimensional array containing the independent variables used. The first 
subscript denotes the independent variable and the second subscript denotes the 
data set. 

YO is an array containing the dependent variables which correspond to the 
independent variables. 

LN is the number of data sets. 

N is the number of independent variables used. 

IP is the degree of the polynomial curve fit. 

These arguments must be assigned a value in a driver program. 

The regression subroutine incorporates several special routines: RDWT, IOWR. 
INVRT, and DMATML. 

RDWT is a general FORTRAN I/O package which is capable of reading and/or 
writing on magnetic tape or drum in FORTRAN or non-FORTRAN format. 

IOWR is an assembly language I/O routine used by RDWT. 

INVRT is a double precision matrix inversion and simultaneous equations solver. 
The call to this routine is: 

CALL INVRT(A,N,M, DETER) 


where : 


A is the input matrix for inversion or augmented matrix for simultaneous 
equations. 
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N is the order of the coefficient matrix. 


M=0 for inversion only; otherwise, M is the number of constant vectors. 

DETER is the determinant of the coefficient matrix. 

DMATML is a double precision matrix multiplication routine with transpose 
options. The call to this routine is: 

CALL DMATML(C,A,B,M,N,K) 


where : 


C is the product of matrix A times matrix B. 

A is the first input matrix, 

B is the second input matrix, 

M is the number of rows of A, 

N is the number of rows of B, and 
K is the number of columns of B. 

Certain variables have dimension changes that vary depending on the number of 
data points (LN), the number of independent variables (N), and the degree of the 
polynomial curve fit (IP). These will be dimensioned as follows: 

DSN(NCL2), DSVHATfNCL, 1 ), DB(NCL,1), X(LN), Z(N,LN), ZBAR(N,IP), 
ZZBAR(NN,IP), SHAT(NCL,NCL), YO(LN), SVHAT(NCL,1), YC(LN), BZ(LN), 
RESID(LN), B(N,IP), A(NN,IP) 


where : 


LN,N,IP are as defined above and 
NCL is computed as follows: 

NCL = N*IP+ N!/[2*(N-2)!] 

NCL2 is NCL raised to the second power. 

NN is N minus 1. 

SMALLZ, SMALZ1, and SMALZ2 should be dimensioned greater than the 
number of data points and the dimension is set in a PARAMETER statement, for 
example, 
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PARAMETER IDIM = 100 


B. Program 

In order to run this program on the UNIVAC 1 108, the deck is set up as follows: 

@RUN,//T JBNAME,320590,UDARBYBIN406,01 . 1 00 
@ASG,T 10.F/1/POS/2 
@FOR,IS MAIN, MAIN 

(Source deck for inputting data) 

@FOR,IS REGRES, REGRES 
(Source deck) 

@FOR,lS RDWT, RDWT 
(Source deck) 

@ASM,IS IOWR, IOWR 
(Source deck) 

@FOR,IS 1NVRT, INVRT 
(Source deck) 

@FOR,IS DMATML,DMATML 
(Source deck) 

@MAP,I AA, AA 

LIB SYS$*MSFC$. 

@XQT AA 
(Input data) 

@FIN 

as illustrated in Figure 1 . 

A complete program listing for the case with linear cross products is shown in 
Table 5 and a flow chart indicating the operation of this computer program is indicated 
in Figure 4., 


C. Output Data 


The desired coefficients of the model equation are outputted in E notation where 
BHAT is the computed intercept of the. fitted polynomial expression and the B 

coefficients are printed out in ascending order of degree (P) with the first P coefficients 
indicating the b coefficients of the First independent variable, the second set of P 
coefficients indicating the b coefficients of the second independent variable, etc. The A 
coefficient designation represents the C coefficients of the various cross products in the 
following order: 
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c ]2 , c 13 , ... , c 1N 

C 23’ C 24’ - > C 2N 
C N-1 ,N 


Also included in the output are the number of independent variables N, the 
degree of the model equation P, and the number of input data points (sets) n denoted in 
the output as LN. Further, the significance of the estimated regression equation is 
indicated by the standard deviation, the multiple correlation coefficient, the F ratio, the 
maximum percent error, and the average percent error. ! 

The values of the dependent variable are also shown for the input data and for 
the computed values obtained from the fitted expression for each set of data input, as 
well as for the residual difference in the input values and the computed values. 

Plotted results are also obtained which compare the input and the computed 
dependent variables as ordinates to the point interval along the abscissa (which is 
normally one of the physical independent variables, but which can be a unit indication of 
each data set in order of input to the program). The computed points are connected by 
straight lines and are plotted with * symbols. The input points are plotted with + 
symbols and the points are not connected with lines. 


D. Illustrative Problem 

For the purposes of illustrating the use of the cross product program and its 
capabilities, the following arbitrary equation was selected for use: 


Y = 9 - 6Xj + 7X 2 + 5Xj 3 

- 4X 2 + 1X 2 2 + 7X ? 3 

- 6X 3 + 7X 3 2 - 6X 3 3 

+ 5XjX 2 + 5X { X 3 - 6X 2 X 3 


(38) 
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TABLE 5. PROGRAM LISTING (WITH CROSS PRODUCTS) 
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TABLE 5. (Continued) 
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c c t: l 

I 

C G i 66 j 

K a 0 S 

coco 

i 

UL1 7 lb 

N Z 

L C 0 0 

w 

001730 

PLRA VG 

OCOU 

R 

001732 

HEfiMAX 

PL' CO 


UfMSbO 

r M r x 

CO DC 

w 

OC i 7 P 5 

k 

C 0 ? rj 

h 

If i b 2 3 

K t S I 0 

C 0 c 0 

rt 

001736 

KHU 

cncu 

H 

001674 

H R 

0 o o o 

h 

0 w 1 3 3 l 

5^ AT 

C o fj J 

.m 

i r r 6 s s 

b^ ALL Z 

r c n r- 

K 

l ' 1 :: 2 i 

S.'-AL? 1 

CCgC 

ti. 

C Q l 1 ft b 

b « A L Z 2 

GDOEJ 

K 

0 C 1 ft 7 b 

SHEG 

loll 

«v 

0 o L 6 7 * 

S •>» L S 

G 0 o L 

H 

t V 1 7 r 3 

S S 

(. CCC 

f. 

l c l 7 3 b 

b 7 ADE V 

r 0 a 0 

a 

C C i / 2 7 

SUM 

00 0 0 

A 

001 sb 1 

SVHAt 

CITE 

K 

1 L l 7 on 

my 

C UUL 

K 

0 C‘ 1 * 7 7 

V A K 

C.pf 0 


L ■. Pft 1 ? 

X 

c 0 0 c 

n 

ti 0 1 6 6 * 

XL 

CLOU 

H 

0 P 1 ft ft S 

XR 



TABLE 5- (Continued) 


DCOO S' 01-1 666 Y B 
C1C00 (t CU0636 ZMR 


CCOO' K _ oor*73 V HAN 

0000 n C006H ? ZZBAR 


‘00ob ’k’ Vo f 5 65 tc 


COoO'R 001667 Yt“ 


0000 R 001702^ ft 


CCILil 
oo uii 
GO 1U* 
GDI US 

00 I :j6 
' C 01 Cl 7 

001 l 0 
GO 1 1 l 
or; I 12 
C L' I 12 


3* 

<1 V 

S« 

6* 

7* 

6 t 

9 ♦ 

10* 


S U b R G U T I n E fttGKES(Z,YO»Lh.N,lP) 

PARAXETEK JOlh«lOC - - 

OOUbLE PRECISION D E 1 1 R 

OCUbLE PRECISION D SN U > . 0 S VH fT <1 2 * t I . PB ( l 2 f 
DIMENSION AOAKlU?! 
r< I - t On PXa* <1 S ) 

DIMENSION FIDM 1 2 ) , FtCt < 1 2 > 

OlMtNSlCN XUS) - 

cl I On Z(3| IS) ,ZeA*(3|3) »ZZflAR(2'3l »5HALLZM0("> • 
•SKaCZI C iCrLMI .SKAL72I HUH) .ShAT I 12, 12» ,Y0« 15) *SVHAH 12,11 . • 


0\ 

-J 


0C1 12 

I 1 * 

• 

YCI 1 S 7 . SZ 1 IS) , 

CO 2 13 

1 2 « 


DATA 1 ZEKO/Q/ 

Cl. 1 1 b 

13* 


data flOy/72h 

c :n f j 

M* 

• 


Go t l 7 

1 S • 


DATA F10X/72H 

CLuT 

1 6* 

* 


0 J 1 

I r * 


DATA < A U A, R Y 1 1 ) 

r:; 143 

1 <* * 


CAUL JOEnT ( 9 * A £ 

C ■ : ; 4 R 

2 V * 


h •. 0 S = I 0 L H 

0 ..US 

i U* 


7. L a C ♦ 

:: r s z 6 

2 1 * 


AH = 1 6 . 

C U ! 47 

22* 


Yfc = 1 . 

DO J 3 0 

2 3* 


Y T = 1 7 , 

c : j 3 1 

2S* 


a R l 1 E ( 6 | 1 0000 J 

C c i j 6 

2 5- 

ItiOOO 

FO- HAT { 1 H 2 * 2bN- 

s. L 1 J 0 

4 v » 

L 

1 M 1 1 A L I L t VAR 

COI 37 

2l* 


M« = N- i 

cr; j <♦ a 

2 if* 


M P = N* JP 

com 

29 * 


DU i !i I *» 1 * K 

C I) J ■« H 

3D* 


U 0 ib -J* 1 » 2 P ... 

C 15 1 H 7 

3 1 • 

1 5 

i ,j ) = o«a 

Cl' 

32* 


DU It I « l • N M 

t)o iso 

J J • 


L> 0 16 J- 1 ► 1 P 

CCi 6Ci 

2<t* 

U 

zzy ak c t * j ) co *o 

COJ 63 

3b* 


0 G 1 I » 1 , L 0 I H 

r L' : 6 6 

it • 

7 

S*Al-LZ ( I >*U.{J , 

CO! / 0 

37* 


DU B i n 1 , r* i p 

CC i 7 j 

i(--« 


DO fa J= 1 ,N I P 

r-,; i j A- 

39» 

6 

Sr. A T ( i , J ) =D . D 

r r.: 

sO • 


i/U 9 I r I » N 2 P 


H J - 

V 

SVHA7UiU-Q*0 

? i’ 0 6 

*! 2 ♦ 

„ 

DO S J-l.LN 

COil 1 

H 3 • 

c 

bZ f 1 1=0-0 

C 0 i 1 3 * 

4 L ♦ 


Y B * K s(5*C 

C 0 2 : *; 

5S » 


Rrt ru ,i. 

002 15 

Hi* 


= 0 

CL'i 1 6 

4 7 * 


SRtSfcC.Q 

G'11217 

Yfc ♦ 


VAK'O.U . 

Ct\ ! ?C* 

** S 1 * 


s y v - o * r 

l;f,, i l 

5 !.’ - 


h - D • 0 


i Values 

point intervals 


U 1 ' 


f 6 H C DP T , 6HF LO 



CT\ 

00 


TABLE 5. (Continued) 


00222 

5 1* 



V Y = 0 « 0 

D0223 

52 * 



55 = c *a , 

CQ223 

53* 

c 


COMPUTE ZflAR 

0022* 

S*« 



00 J 25 I»i ,14 - . ■ ■ ... . . - • • - > 

CC227 

bS. 



00 lib Jo 1 , JP 

CC232 

6 6 * ' 



UO JOa K*l*lt* 

CU2 35 

57 4 


100 

2hAkM,J»= Zb*MI,jU 2<1»K>**J 

CC 23 7 

5C* 


125 

/b^H.Jf 3 ZbAK<!,Jl/LN ’ 

00237 

5 V* 

c 


COMPUTE 2 Z U A H 

CC2*2 

t 06 



K ** 1 ■ 

ct-zia 

6 l • 



00 J 50 1 « I ,NN 

C c 2 4 6 

62 • 

” 


Kt>| f| •■— ■ - - — — 

C32*7 

63* 



HO 150 J»K,N 

C0232 

OH* 



t-G i 30 L- J *LN 

CU2a5 

65 ♦ 


1 3 C 

2ZbAM(i B j)pZZ8AR(rlJ)+ZlI|L)*2(J»L) 

01*2 b7 

66* 


150 

2 2bA« ([ »3 1 eZibAk ( 1 t J 1/LN 

CC257 

67 * 

c 


COhHUTE SMALL Z » S 

OCi* 2 

6*3 ♦ 



i o r t a j o 

C 3 2 *i 3 

6 V * 



i«5lC T e \ N H 05 + 27 l / 2 F 

0 0 2 f»* 

;u« 



I5CCT a 1 

CS**.£, 

7 1 * 



L>0 20C J » 1 « ft 

tl,70 

72* 



CO 2C0 J-l i IP 

:;LV7 j 

7 3* 



LG 175 It - 1 , L h 

C027A 

?*• 


l 7b 

SKALLZ(KJ«2< 1 iSl* 20 AK l I r J ) 

C r. 3 L li 

75* 



CALL lu ; ?GV.K<lilUNIT«]SLCT,SNALLZ t Nf.D5,U*«) 

C::.*o : 

76* 


213 C 

1 bhCT = l5ECT+f,SECT 

Oo jo* 

77* 



KK s 1 

C u j u S 

7 L ♦ 



00 230 i a 1 |NN 

0 r > ic 

7 7* 



KK s* H* ) 

cl j i l 

*!?• 



VQ 233 J » MC . N .... .... ^ ... .... . 

C C J 1 * 

6 1 • 



00 2 2 5 K*l t LI* 

C 7 3 ! 7 

«2* 


225 

SfULLZiK } *Z l I »K 1 *Z1 J ,H >-2ZbA« (j ♦ J) ■ 

0032 1 

fa J * 



CALL ICRpufH 1 ,niNlT«lSECT,SMAULZ,Nfk05,lti«K) 

C0J22 

UH * 


230 

IbtCTa iStCT+NSLCT 

0 0 3 2 2 

as* 

c 


COMPUTE N u fl b L « CF C 0 L U H W S 

0 Vi2*i 

6 6 * 



K h A C T » 1 .... .... _ 

GO 3 26 

a 7* 



00 25Q I«1 iN 

00 33 1 

6»3* 


2SC 

i,mct-nfact*i . . 

0 0333 

BV. 



N F A C T 2 □ 1 

CO.U* 

vo* 



Iy2*N-2 

CC J J b 

V i • 



J F ( N 2 . E 0 . 0 J N 2 e l 

GO 337 

V2* 



00 2 7 S 1*1 i N 2 . . _ * 

CO 3 *2 

Vi* 


27b 

mF AC T 2-MFACT2* i 

C 0 J *f * 

V H * 



KCOLaN* 1 P+ ( WFACT/ 1 2*NFACT 2 1 1 


*5 ♦ 

c 


COMPOTE ELLMtWTS OF THE S h A T M A T K J * 

C-Ji'fb 

V fa * 



l 5E C Ts 1 

003 h 6 

v ; • 



KStCTa<N»i05 ♦ 271/2* 

003*7 

v y * 



I K0» « J 

CG3S0 

9V* 



UO 310 J^l.MCOL 

C-J JS3 

It J«> 



CcLL r u K 0 *< * i 2 t J U M I T , J 5 E C T . 5R A L l 1 , N rt 0 S t J t R R ) 

UO Jb* 

1 13 1 * 



1C0L*1 

0035 5 

lu?* 



1 5t C T 1 ■ i 1 

0035 6 

10 3* 



00 JOS Ko l »HCO l 

CL36 l 

i a* • 



CALL !U«0WH(2 , jUtjll i 1SECT1 ,Si4s,L<t2«rtAuS|*Ei{K> 

00 362 

l J 5 * 



HO JJ” i«i,LU 

GCJoS 

1 Li 0 * 


3 CO 

5n A 1 < ipOiV. ICOL 1 *SHAT< IrtOA * 1C0L )**SPAlZl (I)»5MALZ2M] 



TABLE 5. (Continued) 


C C 3 6 / 

i 1/7 * 


jcol* icaL+1 

c a 3 7 cs 

\ oa* 

305 

1 5tU l s J S E. C T l +NSECT 

CC3/2 

J L’V* 


1 H 0 . V » 1 K 0 A ♦ l 

CCJ73 

l 1 0 * 

3 1C 

lSFCr = j S£CT*NS£CT 

COJ73 

Ml* 

c 

COhF’Otc StfhAt 

Z C 3 7 S 

M2* 


DO 3 2 5 f * 1 * Lh .. — .. . - 

CO 1 * JU 

t 13* 

3 25 

Yt> * tO < I ) 


M '4 1 


Yb*K«reAK/LN 

C C *1 *J 3 

l 15* 


i 5tCT=l 

C fc 4 C *» 

1 J 6* 


UO 15Q KdI.SCCl - - • 

Of? hL 7 

M 7* 


CALL lC.H0AM**l0MTMHCM5MAL2l f h^DS,l£RR) 

C C ■* l c 

l ! « ♦ 


0 0 3 i 0 1 * J « L fJ - - . 

CCI 3 

l M* 


Si/H* 1 IK » 1 MSVHAT IK, M*SWAL21M 1 •MC< I J -Y8AR) 

fSMl 

I2G* 

3 3 C 

0 b V H A T (Ki 1 ) «* 5 V h A T { K » | ) -■ 

Cir.l* 

l 2 l * 

3S0 

ISECTa lb£CT*NSeCT 

C L'H2<J 

3 22* 


12*1 - 

'. v •* i 1 

123* 


00 375 Ml *HC0 l 


124 * 


OQ 375 J=1»MCQl . . . ... - - . -- - - * 


1 75* 


US'* M2) =>5HA1 II .J) 

CC-'MO 

\ tt>* 

375 

- — - - 

0M33 

\2f* 


fcKnhCOL 

c • : j « 3 4 

1 2** 


121*12-1 


12V* 

C • 

FJnO ShAT MATRjn INVERSE 

* •j s 3 s 

l JU* 


f|o(! V . ... . - 

r i' : i j ^ 

13 1* 


CAlL [frVKT < |j$M *KR ,M , DE TER 1 

C C H J 6 

132* 

c 

C 0 -4 P 'J T L d AND A COCrriCUMTS 10* THL ASSUMED KOQCL 

t 1 C H .1 7 

1 3 J* 


CALL OMATML(DBpDSh«DSVHAT *NK t NR , 1 ) 

0 3 << S U 

1 3 ‘i * 


J 2=1 . . 

C Li ** 1 1 

1 35 * 


00 400 1*1 iN 

C l’ H H H 

136* 


00 4 u 0 0*1, IP I ... . 

CDH h 7 

137* 


5 < T ,JMDbM 2 , 1) 

? C 'j v 

1 J t 

■' i 0 

1 2 a H * i .. _ 

r. ; ? 4 5 3 

13V* 


W1-I2 

C U - ^ * 

1 *4 0 » ' 


6 « 1 . ... 

D Lj *i ii S 

Ml* 


u 0 H l o j K 1 » '•* H 

0 D “ 6 U 

M2* 


K » K * l . .... 

D C H s I 

1 ■! 3 * 


UO 4 i n JaK t fg 

ryithS 

I 4 4 • 


& <1 . J 1 * 0 5 1 I 2 \ , l 1 

GG'«6S 

1 45* 

H io 

i2i = m*i 

d o / a 

Lift* 


OHAl =Y&AH 

COS 7 I 

l H 7 » 


00 420 Ul t N 

CL' '<7H 

M <j * 


D(J 420 Ja 1 » IP . - 

C C *i 7 ? 

MV* 

4 20 

UMAT«B^AT-OU.J)*ZBASMtJ) 

COiC 2 

Ml 1 * 


A"j .... .... . 

Of/ ju 3 

Ml* 


DO 42b \ = l *NN 

0 i? £ £J 6 

l 52 * 


K*K + 1 

or s j 7 

M3* 


DO 425 JaKtN 

rj :j *» i 2 

I'M * 

425 

aTz«H*T- A ( 1 t J J *ZZP AK l 1 . J 1 _ - - 

c c ,5 l S 

i55* 


KR l T£ { 6 . IDO 0 6 \ BHAT 

0 U 5*0 

15 6* 

l 0006 

FOhM aT 1 1 flO i Shklrt AT«£ 1 6 * 9 ) . . .... 

0 9S 21 

1 S7* 


no 4 5 n ! * M L h 

0052^ 

Md* 


DO 4b0 J- » »N ... - 

C0527 

Ibv* 


UO 4S0 K«| , IP 

CCS 3 2 

1 60* 

H5C 

H/ | n=h/ M 1+t? \ j. a K J *Z l J , I 1**.R . ... . — 

005 36 

1 6 1 « 


DO 4 7 tj I * 1 , LW 

CCS HI 

l 6 1 * 


K\5 l 


On 

NO 



-4 

o 


TABLE 5. (Continued) 


C0i4 2 

1 6 3 • 

CJ054S 

U 'i - 

CO 5** 6 

1 

00551 

1 6 6 » 

0CS51 

16/* 

0 0555 

l i, R * 

CO 560 

1 6 V * 

C G b 6 1 

1 1 1; • 

0056 2 

1 71 » 

00563 

1 / 2 * 

G C 5 *> 5 

1 / 3 * 

C b S 6 6 

1 7 R ♦ 

CC j / i 

1 75- 

C05/2 

] /* • 

OEy 7M 

1/7* 

CUS/S 

1 /ft* 

CCS/n 

IV 9* 

OOnOl 

l HO* 

00'iC*4 

i S l ■ 

lifJbuS 

1 0 2 • 

P-3607 

1 r.» 3 * 

GC6 i 0 

liM* 

DO* 1 3 

1 d b • 

COM4 

1 a 6 * 

OOo 1 5 

l J 7 * 

006 i 5 

I 

006 l 7 

J iiV • 

C 0 h ?. 0 

190* 

0062 3 

1 9 1 * 

CDc2t 

J V2 * 

C0626. 

1 73 * 

r r . v ■; 

1 i H , 

CUfcJl 

19b* 

C06 J2 

I y 6 - 

C )t>32 

l */* 

CD* J3 

l V4 * 

C 'j 3 6 

i 9V 

ODt-37 

2 0 U * 

CO*“iO 

23 1 ♦ 

C 0 6 'i 1 

202* 

O'JVt j 

2 U 3 * 


2 <j R * 

C 0 6 * 5 

2 j 5 * 

: 0 6 b 2 

2 0* « 

I Z *> :j 2 

2..I / • 

CO *b 3 

C‘J-> • 

0 0 6 b j 

237- 

00*56 

2 i ■'? * 

00*6 7 

2 1 1 • 

*■16 70 

2 l 2* 

C06/2 

2 1 3 » 

C06/3 

2H» 

006 7 4 

2ib - 

"067 7 

216- 

OO/oO 

2 1 / » 

ou ;*ji 

2 16* 


60475 

KK o « K * l - - 

(JJ 475 LiKK.N 

475 UZ ( t > *E*Z I 1 I ♦ A t J , L t • L I J , I ) *Z ( L , I J ... 

c commute: the t points 

OU 1>00 I a 1 * L N - - 

rc { i > oeiha r tyz u > 

i » = yo { i >-yc ( [ > 

««»««♦«£$ 1 CM 1 >•• 2 

5 DC COnT E»uE ..... ------ 

SOM a 0 • 

17b SO l L - J ,Ui • - - - 

rMAA(L»-A 3 5 (<Yb<L»-K<U)/YClU> 

SClSU^SOh + PPAXlL) ------ 

PEMAVGk I SOM/LN > * \ 00 . 

AMaXbPUAX J J ) - 

00 be 5 L=2,LN 

[F|Ai1M-Pf**xUll 506.5C&.50S 

506 A»AA cPl'AA { L > 

SOS C \J N T 1 NUt ■ .... 

PLK^AXaAMAA* 1 0 g • 

1 J 0 b 0 7 L c I 1 L N 

!t«k* A* St kESJC <L y 

507 COrMlNoE / 

C CCkPoTE STANDARD DEVIATION • - • 

,N 2 *C 

00 b l 0 1 ■ 1 « NM - - - - 

rj I sN- 1 

K <1 a T » * f * I p +- t i 2 

GO TO SJ3 - - - -- 

bTA 0 t 7 -SjSTTKK/(LN-KwAT-l)) 

5 J 3 C DjJ J f HvK . - 

C COYOTE THE «UbTK*L.e CORRELATION COE f F J C 1 £f* * * AND FRaTIO 

1 ) 0 5 2 5 ! * I . ♦ U '* ......... 

5f<kSvSF}fcG + (rc<I’-Vi}Aif)*'Z 

5YY»5 YY* i YO W l-TB/iR » **2 ...... 

S5B5S + lYC(rT-YjArt|*IrOU»-rC(m 

5255 rtLb» 5 N£S*kF;Sl 0 Ml *"2 _ _ _ 

Hn A f -SREG/SYY 

Fi<*T JOw(SrfEG/ UHAT-I } »/ISS£S/UN-|C*tATn 

AklTSHA.lOOt'VI STAOi:V t KHAT,FHATtO 

lOOOV F CJ * M A T ( LH0.1VH$TANDAKO‘a£>/IAriONPEj6,9 l HX I 33rtrtUL.TlPLE CORRELATION 
• C 0 E F F F C l E N T = 1 1 6 » V „ H * , H HF KATlO*El 6 .'?j 

AR[Tt'l6»)091U> 

10010 FOkMA T I IMO i 14KB COEFriClENfSj 

akiTEI 6,100111 (Itf r I . J J p J «1 « IFJ 1 ,N> ; 

I 001 1 Format i lh ,St2t .9/ } 

% 7 rtt r e 1 a » j 0 0 1 2 > 

10012 F0» H AT ( iHO, IHHa COErF I C I LNT 5 1 

R = l 

00 b 3 C l-l , NKi 

« = *♦ i .... . . 

DO S30 J=X,N 

'*.K 1 T t I k , 1.00 1 4) A l I ,0 J 



TABLE 5. (Continued) 


0 0 7 0 6 

2 1 9 ♦ 

1 OQ l H 

FORMAT I' IH *£21 ,9 I" 

00707 

220 * 

530 

CONTIMUE 

00 7 i 2 

22M 


AKlTtUi 10 3 2 6 ) 

007 14 

2 2 2 « 

10Q2A 

FOfiMAT(lH0.12X.OHOSStRvED.l2X,6MCOfiPUTeo,l2^.aHHt;SIOUALJ 

C 0 1 1 5 

223* 


00 12 3 3 l 5 » ! i L N 

0 0 7 2 0 

2 2 4 ♦ . 


1..4ITS <6 ,10027 I . YOU S 1 , YC t 15) ,H£SlOi IS) - — 

0 0 7 2 5 

2 2 5* 

12 3 3 

CGfvT I \0f- 

0 0 7 2 7 

226 * 

I 002 7 

FORM at ( Ih * 3L2 1 . ?/ i • - 

007 30 

22 7* 


Li-"! 

CO 7 3 1 

2 2 * 


>1 

-J 

GC7 32 

229 • 


X ( 1 ) * 1 . 

00/33 

23 U * 


00 64V L 1 » 2 , L W - - -- 

CO/ 36 

2 3 1 * 

6 4 9 

X < L ! )-Xi|_i-l)M* 

c o ; h o 

2 3 2 * 


Call 3U U3L (L 1 ,XL . x'4 , Yb » YT *40 ,fLOX .tloy *-ln , x ro 

r c / '4 i 

23 3* 


CALL 3UU3L<L2,XL,KR,YHir’f,3H l FLDX,FLDY i LM,A,YU> 

G 0 ; 4 2 

ZJH* 


*RlT£[6*lGU33) P£tfMAX,PERAVva .. .. 

OG 7 h 6 

23b * 

1003? 

FORMAT ( IHU.22HKAX I*L:< PERCENT E H R 0 R = £ 1 6 * 9 , A X , 3 2 H A V £ H A G E PERCENT ER 

r r» 

236* 


*HQR*C!6«Y> — - 

u C 7 '* 7 

2 3 7 * 

1234 

- C 0 r* T 1 n u £- 

i ' 0 7 -J y 

2^0* 


CALL end JGl3 . • - 

C 0 / 5 i 

2 3 9 * 


STOP 

0 G 7 h 2 

2 bG * 


End 


£f<>) Of COMPILATION: NO V 01 A&NOSTICS, 



<1 


TABLE 5. (Continued) 


4 F U « * 1 S K w *- T I (9 □ A r 

HSO 009- P t / 1 6-1 3 lb 1 l *0 ) 

SUdRGUT2:*e h{J*T Etf T K Y POInT CCC270 

RLUrPft few TRY POifcT COO? 7 i - - 

.‘•-SITE* CNTfiY POINT C0C370 

CLOSE EMHY HOIKf C"C2^ 

CPE* ENTRY POlfxT f ^ a 3 l 

ioKuak Entry point GPCjAO 

ilJKfcL LSi-'.: CuliiuJ DCC^C^i LMME) CGGi*/, 6LAN* cCh^Cw(^) OCCOOC 


*7. ?E FiNAL rE s <K£SCcS IuLOC^j NAME) 

COG 3 10 VS 

UOOH <NERh 2S 
03Cb S .. y L i 

0206 MCH 
C 2 iJ 7 N i 0 7 i 

0 0 l L. '.STOPS 
li J l i . N t k •' 3 * 


STORAGE ASSIGNMENT (bl_OC<, TYPE, ftflM HE LOCATjON, Na:1E) 


u u f) l 

roc 0 7 1 

iOL 

oooo 

000070 

tcor 

ot?et 

000116 

20L 

OGcl 

CCO 1 **G 

2 1 L 

0001 

DQD2<42 

215 6 

<yyc i 

i hi 

22L 

0001 

000? i 1 

2 -a 

coot 

0-30 22 7 

2&L _ ._ 

0001 

0 C 0 2 3 H 

26 t 

OOO l 

000033 

6L 

OOC l 

L*C C J -1 b 

71 

coo 1 

oDOrs? 

Su 

cooo : 

1 OOOQ61 

1 

0 0 0 0 i 

[ 00UQ6S 

1 ulank 

00 0 0 1 

[ 000000 

IC*D 

r - : 2 c : 

i 

; pt-o 

or oo 

C C C | I 1 

IN-JPS 

CCCJ 1 

i C3J066 

I o-" 1 

C3(jO J 

l OOJG** 1 

JPACKT 

0000 J 

[ 0CMG67 

1 T V. 0 

UGcl 1 

[ C :,, JG6 2 

1 i 

0 0 00 1 

t ucO'ua 

[2 

COCO I 

( 00006** 

. J 

COnQ j 

[ 000056 

<NTAQ 

ooao i 

[ OOOOHQ 

kcoe 


00 1 i:o 1 i C •*>♦*****»*•♦•*♦♦*•««**♦««*«* ^ * ♦♦ • 

c o : u J ^ c glnlhal forth an i/o pack a 3 c. caec / r i i a** it , i ■* 70 

CGlijO 3* c H£ i TL KAS TApF Ok DRUM 1 .V FOsMa'I U H NON - r 0 * T R AN FJ fl rUT* 

OflJC HI C ♦ * * * * * 

Ofjlvl 5* SJtfHoJTtNE ROM 

r -p;,3 6 4 Bl«t\5lON AKKMI . I ) . I C »D I 3 J t 1 F l J ( 2 i ) , <00 £ I A I , 1 P A C < T I 4 1 ♦ A N T V,0 ( 3 > 


r 

. ’ z '♦ 

7 * 

£ i T A 

l l F L 0 t l ) »J*=l »2^>/lHl 

, 1 M ?. . 

1 H J , j rt -i t 1 ii b , 1 M 6 

» 1 4 


C 

■; : t, H 

6 • 

12« 1 I . 

2 r l 2 . 2 H 1 3 ► 2 H l H , 2 h t 5 , 

2H 1 6 * 

2 H L 7 • 2r1 I 3 , 2H 1 V ♦ 

2 rl 2 

0 i 2H2 3 . 2*U2 i 2H23 

c 

0 1 !.j *» 

i • 

22rs?k , 

ZkZ'j ,2 H2* , 2H 2 /, * H2B , 

2H29 / 




c 

0 t 'J 6 

l u * 

DATA 

iko3E 1 1 > « i « i * 6 i /ooo * j jonno aoc . o j r !2-jo«.jooo 

30 0 

,000 1 soooooooo * 

C 

? i o 6 

l ) » 

l 

□ r»rj i 

nuOOC 

J ,0 0'? 1 lUOCOV 

00 J 

i ouDM onoooOGoc/ . 

0 

n i l u 

i 2 • 

data 

J |/262 I'lH/ , J 2/ l ?73M h2*l/ 

i 1 3/00 1'Jfl J 1 / 



'j 

l l 1 H 

1 3 • 

j A T a 

1 kLA<«K/6rl 





0 

Oil 1 * 

1*+* 







« 

0 116 

1 b « 

L * r k Y 

k £ 0 T P •* l 1 UMT,kjjl»1 

r ;? « , n 

« » N :■ o 5 , A * '* A i ) 



c 

0120 

I n * 

lOP^M'lOE -7 

^ 1 t k f. 

A» FT’t FJk).AT, 

= 2 . 

ktAi} ;i04-FTN. 

0 

C ! 20 

1 7 * 








G 

012 1 

18* 

'g j T 1 1 

10 






■:• 1 2 2 

l V * 

C..I f •% Y 

k l TCk (1 g K 1 T » u OL . 1 . r '' 

'» 0 5 « A k n A f ) 





TABLE 5. (Continued) 


cm 21 

D0 123 
P 0 1 2 i 
C U I 2 7 
C U 14 0 

001 Jj 
C U 1 3 1 
0JU3 
CD 1 S'* 
<;r i 33 
C C I 3 4 
C L' j J 7 
C 1 J / 

0 0 ! -t0 

: r i '♦ i 
: M s 
01M HS 
CD l 4<i 
00 1 
C C i 5 0 
DO 132 
0013 9 
C 3 13 4 
00 1 3 3 

c : 1 3 s 

3 0 j 3 7 
DGi 6? 
00 j 6 t 
DO 162 
0 .? j ft 3 
CD|q4 
0-W6H 
00163 
CJUb 

cu i67 

30170 
COi fO 
oy 1 7 1 
O 0 i 7 2 
0 0174 
001/6 
00177 
0 0 2 U i 
r ; 4 j 3 
c ■; j j m 
0 0.' 15 

uoz HI 
0021 1 
0 0 2 U 

002 1 3 
002 2 I 
0027 1 
03222 
CJ223 


20 * 

2 1 * 

22* 

2 3 « 

2 + • 

23 * 

26 * 

27* 

2d * 

2 '♦ 

JiJ* 
j 1 « 

J 2 * 

J J* 

11* 

J j * 

}7* 

Hi* 

3V« 

*to* 

1 I * 

•i l * 

4 J • 

4 t • 

4 b • 

'16* 

H 7 - 

4 j • 

4 >■ 

«J| A jhOST tC« 
b :* 
b I • 

32* 

5 3* 

«01 AuNOST l C o 


i aps*oo=>2 ' ft 

<NT*0 ( 1 ) *NAQS » J 1 + t 3 G 
K VT -«0 ( ? > a It! LAN* 9 

K M >'• 0 C J ) • < H T it 0 l I > 

SC TU 10 

* 

ENTHY CLOSE l [UMT, JOP 1 
t T * 0 3 1 

GO TO (6.6,7*23) . jJP S 
I T /. C = 2 * 

I 0 P » b ft 

GO TO 1 Q 

♦ 

ENTRY Cr»EN ( Ioiul i JO p ) 

Jr { JOP.GT « 1 ) Rt TURN a 
I OP” 6 ft 

G u TO Ip 


**3».YfUTe kTm fjehaT. ■*I»«kITe KpN-FTN. 

901L0 POWTHAM ^ecoKO-CONT^OL A0«0. 
OUil.RY CHECKSUM ERROR AOKD. 


J0 (> *4, NO ACTPJrt ,Ht rUHN. 

J 0 P a l * 2 * V*«ITE t'JF 4N0 HZtltiD I'jNlT* 
ab, *> K l T £ tND-OF-FILE ■ 


J0P-2.3. NO ACT IQN iKtTURN* 
- 6 • K L .1 1 H U i U N I f « 


cntuy johoar; icme , Iunit , r array , iekr> 


all St t 0 £ li T 1AL -iit[T£. 
a2i SEQUENTIAL READ. 
aLGtNNlNG Drt’JM-SCCTUX 


$4 « 

3 ’J » 

'3 6 * 
S7* 
&d • 
b 9 ” 

o 1 * 
6 2" 
6 3 « 
6 *4 • 

«. 3 » 

66 - 

6 f * 
6 8 " 
t» i * 
/O. 

7 1 - 
72- 


I r { I CODE . E« » J J I OP”N 
I F 1 1 COOt .0R . 2 ) l OP 
I P A C K T { © ) = t AUDKS 

10 CONTINUE v 

i I ' A C K T ( ! I “ 1 F L D ( I U M T T I ft 

{ P A C < T l 2 | "] iJLAhm 
JP*C/<T r 3 I aO 

l | i P a C < T ( 'I | »KJuE c I OP | ft 

50 70 I 2 J 1 2 1 .20.2] , 22 , 22 ) . I OP 
?C IC-.Om a 1 " i l ♦L3C i KNU - . J I l J J , ft 
IrtF WtUQLS ARRAY lb u1HCIJS]QnLD* 
1C"7(?I«'^!1SMH!.3 CIA«i:A?; ft 

1C.»0(3J- 2 * U-4L0C \ KtiTAOl 2 » 1 ft 

IP aCJCT <ft> « 3*lHl.JCUC.*0(lP ft 

GO Tj 2l 

7h£ Y A R I A 3 L E ArtRAY lb D 1 .11' *5 1 GMCD « 

I P.A C < T | b ) « I 1 +L 0C < A?(RA Y J f» 

lO'.RT ZPACiS T ) ft 

OP.EO.o) RlTJRN ft 

IF ( I 0 P * E • 5 ) Go To l2S.Bj f [T.*0 \i 
ItKRa[PAC<T<N)/l2 ^ 

IKtlLlk.til.Q) bD TO 24 <i 

lMlcRE.0E.3l \iOT0 26 ft 

1 £ E R * 2 ft 

return 

( J = AE Ti.D t t I / 1 1 
IF | ,N00c .E5?. 1 > 50 TO 2b 
tVA* IP AC* I l H IrlPACKt (4 1 t 2.3 9 1 1 
TERrCs] 

RETURN 

n K . T i; I A , 1 0 U > llFACKT(lltl*li6) 

FOR "‘A H ///50n* • • • o*EAECOT 1 go T E i< P 1 N A T E D 
U4i; I/O PACKET J‘j 6U16 I 
STOP 


1U0HESS. 


2) 

22 CALL 
IF t i 


24 


2b 


2 6 
1 0D 


[OMIT IN FI ELD-OAT A. 

OP^KATION CODE TO lOaS. 

FIRST CONTROL .ACCESS. 

SECOND CONTROL ACCESS ^ORL>._ 
THIi^D CONTROL ACCESS *0«l), 
hi«3« AND LUC OF 1ST C7 a ^ORD. 

f» AND LOC OF DATA ^OKDS* 

ft KOOTINE TO ISSUE CXFC «eOu£STi 
ft RETORT aFTca reminding jOntt. 

ft liao*t ,RCTU«ili =2,f?tAlND I unit 

5 Eh^Oe STATUS I s Si PORTION* 
ft UQRNAL OPERATION iNO ERROR* - 
ft TtR*UNATlrni c ft ft U K . 
ft EliO-OF-FlLE OH t NO-OF - T A P E • 

ft N,i FKOtl FTn COnTHOL ^ D R D * 

Nw > HOrl IOmS* 


IN 'RO-'T* 1/0 ROUTINE* // 


ft 



'J 

■f* 


TABLE 5. (Continued) 


tro of compilation; 2 diagnostics* 


i A $ M , IS 2 J •.» i< , 1 0 ;* fi 

a s 1 1 f 

i - 

0 1 COG Cud 27 .00 IH 12 

3 « U C G 0 3 0 i 7 2 11 00 0 0 

** COC0O2 7 H CM 00 U 

5 • 

i i; a 5 

CNO ASM F<^'<OrtS ; NQN£ 




SCll 

aXKS 


0 

.030000 

j 0 A K • _ 

- 1 

.. AO.D,All 

0 

1 0 D 0 0 0 



lO.ri 

0 

CCn Q o 2 


J 

EMU 

2 » X 1 i 



TABLE 5. (Continued) 


G F 0 K * I S I N V R T , 1 U V K T 

USD 0J9-01/ u- I 3 : bt ( • C > ■ - 

SUBROUTINE IMVHT ENTRY P0!M 00336** 


STOrtftGfc USE.0: COO £ ( 1 ) DOO^lZ. OATAlOI CC02u2; BLANK c 0 H H 0 >* ( 2 J OQOOJO 


EXTERNAL REFERENCE* ULOC** NAKEj 


UJU3 NtrtKJi 


STORAGE ASSlGNMt^T (BLOC*. TYPE, RELATIVE LOCATION,. NAME } 


DO Q 1 

oooo I 6 

n . .. 

oao i 

. OCaiHO 

I I *L 

0001 

000 J 

000 I 1 0 


. aaoi 

CC3 162 

t 7 6 (j 

0001 

OOOi 

J ■} 02 6 2 

73 U 

QUO I 

000332 

2 J5 L 

oooi 

Uui3l 

OOQiJ 7 J 

7*L 

0000 

0 OC0002 

A P A X 

0 DO 0 

QdOli i 

[ 0000 H 2 

JNJcX 

OUQO 

L 0 C G l H 6 

I NO 2 

0020 

0000 

1 0GC136 

J 

oooo 

i cyn’+o 

K 

00 DO 


01JQ02 l. 

.1 l 7i 

... oooi . 

,.„„Q00q3H 

1 23Q 

0001. 

000QH2 

1.3.1 G. 

0002b* 4 

lhil 

oooi 

000202 

20Hg 

OOOI 

00022H 

2 152 

U 0 0 3 Q 6 

2H25 

ooo i 

a 0 0 3 H 2 

2 b 0 L 

oaci 

OOOQ7 1 

7*1 

000 1 Hi 

I 

0000 

J OOOLHH 

1C 

OOOD 

I 00Q1H2 

IHD 

oooish 

UJPS 

fJOOO 

j oouuo*< 

I p I V 

QQOO 

l 0C01HJ 

l* 

U'JO its 

L 

oooo 

1 0 0 0 l 3 7 

NN 

0000 

d ooacoo 

SUN 


0 0 1 0 1 

1 * 

00 J03 

i* 

00 luJ 

3 ♦ 

a :• i o j 

H* 

DlsiuJ 

9 • 

0 0 10 3 

6 * 

C 0 10 3 

7 ♦ 

CO 103 

0 ♦ 

CO I OH 

V# 

00105 

JO. 

00 i J 6 

1 !♦ 

0 U 1 d 7 

i 2 * 

Oul 1 0 

1 3 • 

COU 3 

l H * 

0 0 i i a 

IS* 

001 1 6 

1 6 • 

OG l 2 1 

1 7 • 

0 012 2 

i 6 ■ 

0J12S 

19- 

00 1 3C 

20* 

C 0 J 3 3 

2 1 • 

CD i jo 

22* 

G a l 3 7 

2 3« 

001 H 2 

JH • 

COM3 

25* 

00 1 HH 

2o • 

0 'J- 1 M 5 

27* 

0 U I H 7 

2rt. 

C 0 1 5 1 

2 9* 


SUBROUTINE iNVrtTf A ,'I.H.DETEft I 
PARA IE TER 1011**30 

C MATRIX INVERSION A>iO SIMJL f ANEO^S EQUATIONS SOLVER.. ... 

L 4* . v^UT MATRIX f OR I,m*UON OR AUgOmENTLO MATRIX FOR 5 C «E • EQS* U5Q003D 

c K-OROEH of coefficient matrix usooohq 

C "ISO FOR INVER ISON only 

c M»|vOMOERGFCOWSTanTvECTORS -•-“••■■ 

c iaET£r»:)ETLHM 1 NanT ftF caEKFiCIENr MATRIX 

... OOuaLE PRECISluN A I lJ . DETgK >S I GH * AMAA — 

0 I M t NS J On I P l v ( ] 0 J M ) i l NO £ X 1 1 0 I M , 2 ) 

SUV* j ,000 

00 2C J= 1 .N . • .. 

20 I P I V l J J -U 

00 l tJ2 K-l ,N 

A M A X =* 0 . 0 J 0 ... . 

HO 00 76 1 fl 1 »N 

IF I IP 1 V ( I )-I) SD 1 76 »bO . . ...... 

SO 00 7 S J = I i N 

IF U P J V t J ) - i Ub » 73 t 250 _. 

55 I UO * ( J- 11 **♦ I 

IF! ArtAX-OABSiAl [«om 60,75>75 

h Q l Ks i 

1 c = J 

AKaAsOAUSIAIINOI i 

76 CONTINUE 

lJ'|tfUC>«. LPJMlO + 1 



*-0 

Os 


TABLE 5. (Continued) 


0 015 2 

30* 


IF ( I^-IC>90* 1 

00 15b 

3 1 « 

90 


00156 

32 t 


DO 1 l 0 L= 1 » NN 

00 16 1 

3 3- 


I n o * i l - i ) • n ♦ i a 

C 0 1 6 2 

3 M * 


Uq2*(L-1 ) ♦NMC 

C 0 .16 0 

J 3 » 


A r. A A * A { I N 0 1 

c: 1 6*1 

3 6 « 


M I N LJ J s A { I 14 0 2 1 

C 0 1 6 5 

37* 

I 10 

A(1ND2)=AMAX 

CGI 67 

3d o 

1 IS 

i M D t‘ X U » 1 ) « [ K 

cci ?c 

3?« 


IN,0 0AlK*2J-IC 

0 0 1 7 l 

M rj * 


1 1, l- - 1 ; c - 1 ) ♦ * + 1 c 

C ■: 1 7 2 

m i •* 


hhrt 1 ' l *.***/ ■■ - — ■ 

C C t 7 3 

M 2 * 


ULTc-tf = C ETCH* AGAX 

CO 3 7H 

M3* 

1 HU 

A 1 I N 0 ) c 1 • U 0 0 

0 0 17 5 

MM* 


0 J 15 0 L - l i N N 

o : i o c 

nr. • 


1 N D * ( L - 1 ? * M ♦ 1 C 

CCi’Ol 

M 6 * 

ISC 

A u ‘ 1 » A ( I N o 1 / A K A X 

CC2u3 . 

H 7 - 


DO J 5 1 L« i - — 

C C 4 0 6 

M « 


IF 1 L- F C > 1 6 S > l B 1 . 1 i5 

0 » / i l 

MV- 

165 

iND-llC“l)<M + L ( . 

fj ? 1 7 

5 6 • 


A * a A » A ( 1 t\ D ) 

OL 7 1 3 

b 1 * 


A \ I s o . u 0 0 

0 CV 1 4 

S 2 » 


0 0 l H c I * l 1 N n 

✓ ) 7 

5 3* 


1 «\L 31 ( > . - 1 ) • N + L „ . . 

0 f, /tO 

r _> M- 


I rj f; ?. r ( | - 1 > ♦ H + i C 

C :j ? 2 1 

5 5 « 


AMN0)aAlINC-)-A(U02l*AMAX 

C 0 2 2 2 

S 6 ♦ 

U'C 

COM IKUE 

C C 2 4 H 

5/- 

1S1 

continue 

C 0 4 2 6 

bb , 7 

1*2 

COM I r v vi 0 

t- r 2 :, :: 

b / * 


L'O 23b 1 = 1 * N 

n 0 j 3 

6 L « 


L »I + i " l 

C M J r 

O i * 


iM-UutAU, 1 1 

GC23S 

A 2 » 


IC*lfcftLAlLt2) 

C L 2 5 6 

c3- 


IK 1 IM-1 C 121D. 235 .? 10 

C'L‘2 ^ 1 

ta M ■* 

2 1 0 

0 0 2 3 0 K = l * N 

CC 2vH 

6 -» - 


I N C- a M K - 1 ) * N + K 

OC2MS 

O 6 * 


1 V|J 2 = ( I C - 11 • U + K 

f: 0 2 n 6 

fa 7 * 


AH A X * A ( l NO > 

re? 1 '/ 

(ib< 


A l 1 NO ) * A I 1K’02 > 

00 2 5 0 

<: 9 ’ 

2 3 C 

A l I r > 0 2 1 a A r A X 

C 0 2 5 2 

/L‘» 

235 

L U M l N U £ 

C02SH 

7 1 * 


IsCTF^sSlGN-DtTtR . - 

0 1 2 5 5 

7 2 - 


£ 1 0 X N 

C c >: £■ 6 

; 3 * 

2 5 C 

... 

c :. 2 6 7 

7 M * 

2S'j 

i< r T 0 K JV 

0 t 7 1 0 

7b* 


t nO 


END CF COMPRArjOH: 


NO QJ AGNOSTICS 



TABLE 5. (Continued) 


i»foR,r 5 ut< Ain-. , umaTml 

HiO CiC9-0i/ 16-13 iSl 1,0) - ... . ~ 

bUiiSCUMNE. D I'. A T M L EnTRY POJnT 0001*0 


S T 0 R f • G E USED; COuEM) DOG 1 76 ; 0 AT A C 0 ) OOOQS3; BLaNK cC««0ni21 000000. 


LATERAL REFERENCES (BLOC*, NAME) 
00U3 NLRR3S 


STOr/ GL ASSIGNMENT (BLOCK* TYPE t RELATIVE LOCATION, NAME), 


GOO 1 

0ur;c3i 

IL 

DOG l _ 

ordj}i i 

I 30 G 

0DM1 

OD0D67 . 1 3&Q 

QO.D L 

a_gpm2.j.H3G 

000) 

DOOQ 3S. 


00 r J 1 

o J 0 0 H 6 

3 L 

OOU 1 

CCOcS ? 

RL 

aoao 

0 0 U U 0 Ci C tu 

ooon 

I COOCO* 1 A 1 

ooao ; 

I OQQQO* 

I A2 

OUOU. 1 

[ Lf IGGO* 

I A 3 

OCQG ; 

I CCOrC 7 

I & L 

qocd 

1 L (j U w I U lb2 

oono 

l 0G00L3 I B 3 

oaao ; 

I 000003 

IH 

oucr: i 

t acooo 1 * * * 

I N 

CGCJC 

occo? * 

I NJPf 

coco 

I OOOD1* LA 

coco 

i ooqdl* te 

aoao i 

; coaou 

LC 

C-OOG 1 

r uocou 

L M 

dduO ; 

1 0?0r?l7 

LN 

oooc 

1 OQOOlR LP 







a l i o i 

I » 

Cl' \ L 1 

2* 

SClul 

3* 

C? o i v» 1 

R % 

0 P J 0 1 

h ♦ 

Ct'ic* 

A * 

t-r i u i 

V* 

COIC 1 

b * 

3C10J 

v ♦ 

6 C 1 U 1 

10- 

0 0 I 0 1 

i l * 

CO 1 LI 

IJ?* 

DC 101 

1 3 * 

CC1 0 l 

IM* 

C C 1 Q l 

1&« 

Ci.' I L \ 

1 E ■* 

L C ! u 1 

1 ! * 

0 0 1 1 3 

1 U ♦ 

C0I51 3 

19* 

0 L 1 0 1 

*U* 

CL-IOR 

i l * 

l r; ! 0 H 

f.i * 

0 U 1 Li * 

23* 

CL’ l US 

2*t* 

CVtU 

■it . 

C(MU7 

Zb* 

OCJI 10 

2 t * 

coin 

23* 

COM3 

2 V* 

CPU 4 * 



SUBROUTINE DKATML (Ci A .& »* »N iK ) 
c* abstract 

c b L N t R A L MATRIX MULTIPLICATION KOUTjNL V- I Th TRANSPOSE OPT^S 

C 6 Fi t K E » M IS The wvJ N.-jE R OF r\G*5> QF IA) -• - 

c N IS The MjMcjER U F NOV'S OF ( ti I 

C K IS ThE MHaE* OF COLUMNS UF <b) OR ' B ) T • - 

c TRANSPOSE OPT IONS /> l COMROLLLD 6Y The SI6NS of M and n. 

c The following products MAY BE OBTftiNEO 

C <r)3(AH6J M A i\ 0 N POSITIVE 

c t c J - f ft J t i tj > . m negative for UiT.~. 

C <C>*rAltG)T N NEGATIVE FOR 1 B > T 

c. (ClMiMUjlT H AND N NEGATIVE 

L V. (-: E ri E T INDICATES TkANSPOSL 

c IF M JS NEGATIVE, h t& THE NUMBER Of K&itS OF IA)T ... 

c IF N IS NEGATIVE, N IS Tut MUhUEK OF RO*S GF I3)T 

C - - 

C* OUTPUT ARGUMENT * C 

dimension c i n - - - 

c 

C ♦ IN PUT ARGUMENTS' A. H » M i N » K 

DIMENSION A l l U 0 U > 



OClhLt PRECISION C B t C • A , 8 


I A 3 ’* i 

Ills I AGS 1M 

l k = 1 A * S I N J 

!F|ti . L 7 « 3 J (i 0 TO . 1 

I A I ■ I M 

i A ? 5 l 


-0 

•-J 



<1 

00 


TABLE 5. (Continued) 


C C J I s 

3 1 * 


GO TO ? " 



C c J I 6 

J 2 * 

1 

1 M * 1 


CC J i 7 

33 * 


I A 2 E I N 


c o l ;; o 

iH * 

2 

I M N * tT 

0)G0 TO 3 

0 0 U 2 

3s * 


181 = 1 


00 J 2 3 

J 6 « 


. 1 0 2 = I N - 

_ „ _ _ , 

00 J 

3 7- 


GO TO H 


OUlia 

33* 

3 

I e ] »K 


CDl/e 

39* 


1132*1 


CC] 27 

HQ* 

H 

DO 7 CM=J 

.IK 

C0 1 32 

H 1 • 


1 1 s L M 


G r ’U3 

H 2 « 


ia:.*i 



C L J j 1 

H J * 


UO 6 L P = 

• K 

C C 1 3 7 



CO * 0*GuO 

CC 1 0 

h b * 


L A = I A 3 


CO] H i 

H6 * 


L ^ c I b 3 


C “ 2 

H / * 


00 s t.N = l 

. IN 

C C 1 ' *5 



CO * CO 4 

A t L A ) * 6 1 1 b ) 

C C 1 6 

4 V' 


LasLa+jAI 

C c I S 7 

bC* 

S 

L3c.u& + :ai 

Co ' “ i 

b 1 * 


C t L c ) a CD 

C c ! S 2 

b < * 


L C * C C * IK 

C 0 U> 3 

5 3* 

6 

I ■& 3 « I 6 3 ♦ 1 b 2 

C C l 3 Li 

bH« 

7 

1A3MA3+IA2 

00 i >7 

bb* 


K L T 0 R M 


C c i * Li 

b fc * 


C fvD 



fcNU OF COMPILM 10*: NO DIAGNOSTICS, 



TABLE 5. (Continued) 


tKAP'f A A, A A 

*AP 1 /16- 1315 1 -< 

I* L I 8 $YSs*tf5FCs. 


aodmss Lihj': 
STAKTih,c, AD0KE5S 

*0F!l.5 DLClhAL 


CC.lDOC 033*60 
C 3 3 5 *4C 

1^601 lBA*K 


OtCOOO C5SC12 
66*7 OB ANK 


SEGMENT ?MlN OPlOOO Q3356C CSOOQO. 0550 12 


NS‘.TCl/FO H 
fcpsillt S/hSFCStS 
f.LXF tS/f SF CS7 
A L r Jvj/ F LRb ] 
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Figure 4. Operation of the computer program (cross product). 
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Figure 4. (Continued) 
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Figure 4 . (Continued) 
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Figure 4. (Continued) 
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Figure 4 . (Continued) 











































From equation (38) the following set of input data (Table 6) was developed, 
containing exact dependent variable values (to three decimal places) for arbitrary values of 
the three independent variables. In as far as was possible the problem selected in Appendix 
A was duplicated for comparison purposes. 


TABLE 6. INPUT DATA, DEPENDENT VARIABLE VALUES 


Data Point Number 

Y 

X 1 

X 2 

X 3 

1 

6.759 

mm 

mm 

0.4 

2 

8.251 

^19 

iiil 

0.2 

3 

4.256 

0.2 

m 

0.6 

4 

4.892 

0.2 

0.6 

0.5 

5 

13.165 

0.9 

0.0 

0.5 

6 

9.669 

0.7 

0.0 

0.6 

7 

6.070 

0.3 

0.9 

0.7 

8 

12.588 

0.9 

0.3 

0.6 

9 

5.749 

0.5 

0.6 

0.7 

10 

8.561 

0.3 

0.8 

0.2 

11 

11.333 

0.2 

0.9 

0.0 

12 

10.589 

0.8 

0.5 

0.6 

13 

16.117 

0.9 

0.8 

0.3 

14 

3.525 

0.1 

0.4 

0.7 

15 

2.333 

0.0 

0.5 

0.8 


The input for this sample program was read in through namelist NAM1 as follows: 


Col. 2 
$NAM1 

Z = 0.5, 0.4, ...,0.8, 

YO = 6.759, 8.251, ..., 2.333, 
LN = 15, 

N = 3, 

IP = 3, 


The computed results are shown in Table 7 for this example problem and the 
plotted results are shown in Figure 5. 
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TABLE 7. OUTPUT DATA LISTING FOR ILLUSTRATIVE PROBLEM (CASE WITH CROSS PRODUCTS) 
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. 9 669C00C3 + 0 1 » 96690002 7 + 0 I ’. 230416579-06 

. 60699999 9 fO l *607000053+0 1 - * 5 3 6491903-06 
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Figure 5. Comparison of input and computed values for illustrative 
problem with cross products. 
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The plot graphs the point intervals on the X axis against the exact and computed 
dependent variables. The exact dependent variables are plotted with (+) and the points 
are not connected by line segments. The computed dependent variables are plotted with 
(*) and the points are connected by straight line segments. 

inspection of the digital and plotted results reveals that the computed regression 
relation yields an almost exact representation of the input data for this arbitrary case. 
However, experience with several sets of physical data have generally yielded poorer 
results evidently due to the larger matricies being manipulated and due to unknowns as 
to whether true linear cross coupling exists in the physical process. 
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